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SUMMARY 


A  literature  survey  was  conducted  to  determine  the  state  of  the  art  of 
measuring  and  predicting  aerodynamic  characteristics  of  oscillating  air¬ 
foils.  Results  of  this  survey  are  presented  in  Appendix  I  as  a  correlation 
and  tabulation  of  airfoil  and  finite  wing  experimental  investigations.  An 
extensive  bibliography  resulting  from  the  literature  survey  is  also  pre¬ 
sented. 

Aerodynamic  forces  on  a  two-dimensional  NACA  0012  airfoil  oscillating 
sinusoidally  in  pitch  were  measured  by  two  techniques.  The  forces  were 
obtained  from  pressure  measurements  and  by  means  of  strain  gage 
balances.  Pressure  measurements  were  made  on  the  airfoil  oscillating 
in  pitch  about  the  quarter -chord  point  at  various  mean  angles  of  attack. 
Strain  gage  balance  readings  were  obtained  for  models  with  pitch  axis 
located  at  25,  37,  and  50  percent  chord  points  oscillating  about  various 
mean  angles.  Direct  force  measurements  were  employed  in  an  effort  to 
obtain  drag  data. 

Test  results  obtained  by  the  two  measuring  techniques  exhibit  excellent 
agreement  over  the  test  range  of  oscillating  frequencies.  At  low  mean 
angles  where  the  instantaneous  angle  of  attack  does  not  exceed  the  steady 
state  stall  angle  of  attack,  the  data  compare  very  well  with  incompressible 
theory.  At  higher  mean  angles,  the  pitch  oscillations  were  found  to  in¬ 
crease  the  stall  angle  of  attack  with  corresponding  increase  in  the  normal 
force  and  pitching  moment  coefficients.  Mean  values  of  drag  were  found 
to  increase  with  increasing  oscillating  frequency.  The  oscillatory  ampli¬ 
tude  of  drag  tended  to  decrease  as  oscillating  frequency  increased. 

Instantaneous  pressure  distributions  are  presented  for  representative 
oscillating  conditions. 
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FOREWORD 


The  results  from  the  oscillating  airfoil  tests  are  presented  in  this  report. 
The  project  was  performed  under  Contract  DAAJ02-67-C-0017  (Task 
1F162204A13903)  under  the  technical  cognizance  of  Patrick  Cancro, 
Project  Engineer,  U.  S.  Army  Aviation  Materiel  Laboratories. 

The  cooperation  and  assistance  of  Stanley  E.  Pearson  of  NASA -Langley 
Research  Center  with  the  pressure  instrumentation  are  gratefully 
acknowledged. 
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INTRODUCTION 


In  the  past,  considerable  effort  has  been  expended,  both  theoretically  and 
experimentally,  to  determine  the  aerodynamic  characteristics  of  oscilla¬ 
ting  airfoils.  A  survey  of  the  literature  in  this  field  was  undertaken  to 
determine  the  state  of  the  art  and  the  extent  of  previous  experimental 
investigations.  The  primary  emphasis  of  the  survey  was  placed  upon 
experimental  investigations  of  two-dimensional  oscillating  airfoils.  An 
attempt  to  correlate  the  results  of  these  investigations  was  undertaken. 

Results  of  the  survey  indicated  two  major  problem  areas  in  experimental 
oscillating  airfoil  investigations:  that  of  obtaining  the  desired  motion  and 
that  of  measuring  the  aerodynamic  forces.  In  regard  to  the  first  problem, 
the  motions  of  primary  interest  in  the  pact  have  been  pure  sinusoidal  rota¬ 
tion  about  the  pitch  axis,  pure  sinusoidal  translation,  and  a  combination 
of  those  two.  One  reason  for  interest  in  sinusoidal  motion  was  that  it 
could  be  compared  with  existing  theory.  Recently  there  has  been  an 
interest  in  a  combined  rotational  and  translational  motion  which  would 
simulate  the  motion  experienced  by  a  section  of  a  helicopter  rotor  blade 
as  it  rotates.  There  were  no  experimental  results  uncovered  during  the 
literature  search  dealing  with  this  particular  problem.  In  the  past,  some 
rather  ingenious  devices  have  been  used  in  an  effort  to  overcome  the 
problem  of  measuring  the  aerodynamic  forces.  The  methods  used  fall 
mainly  into  two  categories:  that  of  measuring  forces  directly  and  that  of 
measuring  pressures  and  obtaining  forces  by  integration.  Both  methods 
have  their  advantages  and  disadvantages.  In  order  to  investigate  measur¬ 
ing  techniques,  an  experimental  test  was  undertaken  using  sinusoidal 
motion. 

Aerodynamic  forces  on  an  NACA  0012  airfoil  oscillating  sinusoidally  in 
pitch  were  measured  by  two  techniques.  The  forces  were  obtained  from 
pressure  measurements  and  by  means  of  strain  gage  balances.  By  using 
two  methods  of  obtaining  the  oscillating  forces,  the  advantages,  disad¬ 
vantages,  and  limitations  of  the  methods  could  be  compared  while  keeping 
the  airfoil  profile,  test  installation,  and  testing  techniques  consistent. 
Pressure  measurements  were  made  on  the  airfoil  oscillating  in  pitch 
about  the  quarter-chord  point  at  various  mean  angles  of  attack.  Strain 
gage  balance  readings  were  obtained  with  the  model  oscillating  about 
pitch-axis  locations  of  25,  37,  and  50  percent  chord  for  various  mean 
angles. 

Results  obtained  from  the  two  methods  are  compared  with  each  other  and 
with  two-dimensional,  incompressible  flow,  oscillating  airfoil  theory. 
Testing  at  different  pitch-axis  locations  allowed  an  extension  of  the  com¬ 
parison  with  theory  and  correlation  with  the  work  of  other  experimenters. 
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Mean  angles  were  chosen  so  that  in  some  cases  the  airfoils  were  oscilla¬ 
ting  in  and  out  of  the  stall  condition.  There  is  no  satisfactory  theory  for 
this  condition,  but  it  is  a  situation  that  arises  in  helicopter  rotors  where, 
at  certain  azimuth  locations,  portions  of  the  rotor  blade  exceed  the  stall 
angle  of  attack.  While  the  Mach  number  and  Reynolds  number  employed 
in  the  experimental  investigation  may  be  somewhat  low  for  application  to 
helicopter  rotor  theory,  the  oscillating  characteristics  presented  for  the 
stall  region  and  the  instantaneous  pressure  distributions  obtained  may 
offer  some  insight  into  the  rotor  problem. 
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LITERATURE  SURVEY 


A  literature  survey  was  conducted  to  determine  the  state  of  the  art  of 
measuring  and  predicting  aerodynamic  characteristics  of  oscillating  wings 
and  airfoils.  During  the  course  of  the  literature  survey,  it  became 
evident  that  considerable  effort  has  gone  into  the  investigation  of  unsteady 
aerodynamics.  A  complete  coverage  of  this  field  is  beyond  the  scope  of 
this  report.  Therefore,  the  primary  thrust  of  the  survey  was  directed 
toward  determining  what  had  been  done  in  the  area  of  two-dimensional 
oscillating  airfoil  investigations.  This  information  was  needed  to  provide 
a  background  for  the  experimental  investigation  of  techniques  for  deter¬ 
mining  the  aerodynamic  forces  on  an  oscillating  airfoil.  The  scope  of  the 
survey  was  wider  than  that  dictated  by  the  ultimate  experimental  effort, 
in  order  to  provide  the  following: 

1.  A  bibliography  covering  the  entire  topic  of  oscillating  wings 
and  airfoils  (subsonic,  compressible,  nonuniform  flows,  etc. ). 

2.  Background  information  on  the  experimental  techniques  of  un¬ 
steady  aerodynamic  measurements  employed  by  various  investi¬ 
gators. 

3.  Theoretical  means  for  predicting  results. 

Results  of  the  literature  survey  are  presented  in  Appendix  I  and  the 
Selected  Bibliography.  Appendix  I  presents  an  attempt  to  correlate  the 
results  of  various  low-speed  investigations  of  oscillating  two-dimensional 
airfoils  and  a  summary  of  airfoil  and  finite  wing  experimental  investiga¬ 
tions.  Due  to  the  large  number  of  parameters  involved  in  the  investigations 
(oscillating  frequency,  oscillating  amplitude,  mean  angle  of  attack,  test 
Reynolds  number,  airfoil  profile  and  pitch-axis  location)  and  the  type  of 
data  reported,  it  was  quite  difficult  to  make  any  direct  comparisons  or 
correlation.  All  reports  resulting  from  the  literature  survey  which  are 
not  included  in  the  summary  of  Appendix  I  are  listed  in  the  Selected  Bibliog¬ 
raphy. 
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EXPERIMENTAL  PROGRAM 


An  experimental  program  was  conducted  to  measure  the  aerodynamic 
forces  on  a  two-dimensional  airfoil  undergoing  forced  sinusoidal  oscilla¬ 
tions  in  a  wind  tunnel.  The  forces  were  determined  by  means  of  pressure 
measurements  and  direct  force  measurements  using  strain  gage  balances. 
A  NACA  0012  profile  was  chosen  for  the  airfoil  so  that  results  could  be 
compared  with  previous  experimental  investigations. 

DESCRIPTION  OF  APPARATUS 


Test  Facility 

Tests  were  conducted*  in  a  two-dimensional  open  circuit  wind  tunnel.  This 
tunnel  has  a  capability  of  110  mph  velocity  through  a  1. 5 -foot  by  3.  875- 
foot  test  section.  These  dimensions  are  one -half  scale  of  the  3-foot  by 
7.75-foot  test  section  of  the  University  of  Maryland  7.75-foot  by  11 -foot 
low-speed  wind  tunnel.  The  velocity  distribution  in  the  test  section  is 
very  uniform,  with  a  variation  of  less  than  0.  5  percent  of  the  mean  velocity 
in  the  test  region.  The  turbulence  factor  as  determined  by  a  4.  5-inch 
sphere  is  1. 08. 

Oscillating  Mechanism 

An  oscillating  drive  mechanism  was  designed  to  meet  the  following  require¬ 
ments: 

1.  Oscillate  the  model  in  pure  rotation  (pitch). 

2.  Oscillate  the  model  in  pure  translation  (heave). 

3.  Oscillate  the  model  in  a  combination  of  rotation  and  translation. 

4.  Be  adaptable  to  the  1.  5-foot  by  3.  875-foot  wind  tunnel  and  the 
3-foot  by  7.  75-foot  test  section  of  the  larger  wind  tunnel. 

The  literature  survey  established  the  following  design  capabilities: 

1.  Rotation 

a.  Frequency  range,  2  to  30  Hertz 

b.  Amplitude  of  10  degrees  from  mean  angle  of  attack 
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2.  Translation 


a.  Frequency  range,  2  to  30  Hertz 

b.  Amplitude  of  5  inches 

These  design  requirements  and  capabilities  were  greater  than  the  dictates 
of  the  immediate  experimental  program.  Details  of  the  design  are  pre¬ 
sented  in  a  separate  report.*  A  brief  summary  and  principles  of  operation 
are  presented  here  for  convenience. 

In  order  to  meet  the  requirements  of  pure  rotation,  oure  translation,  and 
combined  motion,  the  oscillator  was  constructed  with  separate  pitch  and 
heave  shafts  with  the  capability  of  gearing  them  together  for  combined 
motion.  The  shafts  have  flywheels  on  each  end  with  adjustable  crank  pins 
driving  Scotch-yoke  arrangements.  Thus,  the  two-dimensional  model  is 
driven  from  both  ends.  This  eliminates  twist  of  the  model  that  may  exist 
if  it  is  driven  from  one  end  only.  The  flywheels  were  designed  to  main¬ 
tain  a  speed  variation  of  less  than  2  percent  at  the  top  oscillating  speed  of 
30  Hertz.  The  crank  pin  had  a  maximum  travel  of  5.  0  inches  from  the 
centerline  of  the  shaft.  This  travel  was  sufficient  to  provide  the  design 
rotation  and  translation  amplitudes.  The  pitch  and  heave  shafts  were 
driven  by  a  10-horsepower  variable-speed  motor  drive.  This  drive  pro¬ 
vided  oscillating  frequency  ranges  of  2  to  15  Hertz  and  4  to  30  Hertz. 

The  oscillating  mechanism  was  designed  to  be  mounted  under  the  test 
section  of  the  small  two-dimensional  tunnel.  The  forward  or  heave  shaft 
was  directly  under  and  parallel  to  the  model  pitch  axis.  The  heave  motion 
was  transmitted  by  a  vertical  shaft  attached  to  the  Scotch -yoke.  The  rear 
or  pitch  shaft  was  located  18  inches  downstream  from  the  heave  shaft. 

The  pitching  motion  was  transferred  to  the  model  by  means  of  a  vertical  shaft 
connected  to  the  Scotch-yoke  mechanism  and  a  pitch  arm  connecting  the 
vertical  shaft  to  the  model  pitch  axis.  In  order  to  insure  sinusoidal  motion, 
the  pitch  arm  was  fitted  with  a  cam  follower  which  rode  in  a  horizontal 
slot  attached  to  the  top  of  the  vertical  pitch  member.  The  general  arrange¬ 
ment  of  the  drive  mechanism,  1.5-foot  by  3.  875-foot  wind  tunnel,  and  model 
is  shown  in  Figure  1.  Figure  2  illustrates  the  details  of  the  pitch  mecha  - 
nism. 

Model  Construction  and  Instrumentation 


There  are  two  main  considerations  in  the  design  of  airfoils  for  forced 
oscillation  testing.  First,  the  weight  must  be  held  to  a  minimum  to  re¬ 
duce  the  inertia  loads,  and  second,  the  structure  should  be  rigid  to  keep 
deflections  as  small  as  possible.  These  two  requirements  are  in  opposi¬ 
tion  to  each  other,  and  a  reasonable  trade-off  between  the  two  is  required. 
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In  an  effort  to  meet  these  requirements,  it  was  decided  to  construct  the 
models  from  plastic  foam.  For  this  purpose,  a  wooden  airfoil  of  NACA 
0012  profile  having  a  1-foot  chord  and  a  1.  5 -foot  span  was  constructed. 

A  plaster  mold  was  cast  from  this  model.  With  this  mold  it  became  a 
simple  matter  to  make  models  using  a  foam-in-place  rigid  plastic.  The 
result  was  a  lightweight  rigid  model  having  a  very  smooth  surface*  The 
first  models  made  by  this  method  were  somewhat  unstable.  That  is,  the 
profile  shape  changed  after  a  period  of  time.  This  problem  was  overcome 
by  allowing  the  models  to  cure  in  the  mold  for  several  days.  Several 
models  were  cast.  Some  were  used  for  checking  out  the  method  and  im¬ 
proving  upon  the  technique.  Four  were  constructed  for  the  experimental 
investigation.  The  models  used  during  the  tests  consisted  of  one  pressure 
model  and  three  force  models.  These  models  are  discussed  in  the  follow¬ 
ing  sections: 

Pressure  Model 


Model  weight  is  not  as  serious  a  problem  for  the  pressure  model  as  it 
is  for  the  force  models.  Whereas  the  inertia  loads  arising  from  the 
model  weight  do  impose  a  burden  on  the  drive  mechanism,  they  do  not 
affect  the  forces  obtained  from  the  integration  of  the  pressure  readings. 
Therefore,  the  limiting  weight  factor  is  the  capacity  of  the  drive 
system.  This  is  fortunate  since  other  problems,  such  as  mounting 
pressure  transducers,  arise  in  conjunction  with  the  pressure  model 
that  tend  to  increase  the  inertia  loads. 

In  determining  forces  from  pressure  measurements,  ten  orifices  on 
both  upper  and  lower  surfaces  were  considered  to  be  a  minimum  num¬ 
ber  for  the  accuracy  desired.  If  forces  are  the  only  interest,  then 
differential  pressure  transducers  could  be  used  and  only  ten  would  be 
required.  However,  if  instantaneous  pressure  distributions  are 
desired,  twenty  transducers  would  be  needed.  Since  transducers  and 
their  associated  instrumentation  are  expensive  and  such  large  numbers 
are  difficult  to  mount  in  small  models,  it  was  felt  that  any  method  of 
reducing  the  number  of  transducers  required  was  worthy  of  considera¬ 
tion. 

One  method  of  reducing  the  number  of  transducers  required  that 
appeared  promising  was  the  use  of  a  scanning  valve  with  a  single 
transducer.  This  would  require  pressure  tubing  running  from  each 
orifice  on  the  model  to  the  scanning  valve.  To  employ  such  a  method 
would  require  a  knowledge  of  the  pressure  attenuation  and  phase  lag 
associated  with  the  system.  In  order  to  acquire  this  knowledge,  an 
investigation  of  the  effects  of  tubing  on  the  remote  reading  of  oscilla¬ 
ting  pressures  was  undertaken.  The  report  of  this  investigation  is 
presented  in  Appendix  II.  It  was  concluded  from  the  investigation  that 
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it  would  be  impractical  to  determine  the  instantaneous  local  oscilla¬ 
ting  pressure  from  remote  readings  using  this  method. 

For  this  investigation,  it  was  desired  to  obtain  instantaneous  pressure 
distributions.  Twenty-four  differential  pressure  transducers  of  the 
NACA  type  described  by  Patterson2  were  used,  twelve  on  the  upper 
surface  and  twelve  on  the  lower  surface.  A  Statham  Model  PL131TC 
transducer  was  employed  to  obtain  the  leading  edge  pressure,  result¬ 
ing  in  a  total  of  25  pressure  orifices.  The  range  of  these  transducers 
was  ±2  psid.  An  aluminum  block  was  machined  to  receive  the  NACA 
type  transducers  in  two  rows  of  twelve  each,  stacked  one  transducer 
over  the  other,  connected  on  the  reference  pressure  side  to  a  com¬ 
mon  reference  pressure  passage.  The  transducers  were  placed  in 
the  block  and  secured  by  a  cover  plate  which  contained  short  tubes 
leading  to  the  pressure  orifices  on  the  model  surface.  Details  of 
this  installation  are  presented  in  Figure  3.  This  block  containing 
the  pressure  transducers  was  mounted  at  the  midspan  point  of  the 
pitch-axis  shaft.  The  shaft  and  the  transducer  block  were  then 
positioned  in  the  mold,  and  the  model  was  cast  around  them.  The 
tubing  for  the  pressure  orifices  was  then  worked  down  to  the  model 
contour.  Thus,  the  transducers  were  embedded  in  the  model  and 
and  could  not  be  removed  without  destroying  the  model  itself.  The 
pressure  orifices  were  located  at  0,  0.75,  2,  5,  10,  15,  20,  25,  35, 
45,  60,  75,  and  90  percent  chord.  The  transducer  for  the  leading 
edge  orifice  was  mounted  in  the  model  separately  from  the  other 
transducers  but  was  connected  to  the  common  reference  pressure. 

Prior  to  installing  the  transducers  in  the  model,  the  system  was 
checked  for  leaks  and  feedback  through  the  reference  pressure 
passage.  That  is,  the  system  was  checked  to  insure  that  oscillating 
pressures  at  one  transducer  did  not  interfere  with  pressures  recorded 
on  adjacent  transducers.  This  was  checked  by  connecting  all  orifices 
except  one  to  an  oscillating  pressure  (motor -driven  piston)  of  the 
maximum  anticipated  test  pressure  and  observing  the  output  signal 
of  the  remaining  transducer.  The.'e  was  no  discernible  interaction 
observed. 

Commercially  available  carrier  equipment  was  used  with  the  trans¬ 
ducers  to  drive  recording  galvanometers.  Two  recorders  were  used, 
one  for  the  upper  surface  pressures  and  the  other  for  the  lower  sur¬ 
face  and  leading  edge  pressures.  The  frequency  response  of  the 
galvanometers  was  flat  up  to  60  Hertz.  TTie  frequency  response  of 
the  transducers  was  equal  to  or  betier  than  the  galvanometers. 
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Force  Model 


Force  models  were  constructed  with  pitch-axis  locations  of  ?5,  37, 
and  50  percent  chord.  Construction  of  these  models  was  relatively 
simple.  A  spar  was  machined  from  1/2 -inch  aluminum  plate, 
located  in  the  mold  to  provide  proper  pitch-axis  location,  and  the 
plastic  foam  was  poured  around  the  spar  to  form  the  model.  Each 
spar  was  designed  so  that  the  combined  weight  of  the  spar  and  plas¬ 
tic  foam  would  locate  the  model  center  of  gravity  very  nearly  on  the 
pitch  axis. 

To  measure  the  aerodynamic  forces  on  these  models,  two  strain 
gage  balances  were  designed  and  constructed.  The  balances  were 
mounted  one  on  each  end  of  the  model.  Each  balance  measured  the 
model  normal  force,  axial  force,  and  pitching  moment.  These 
measurements  were  accomplished  by  means  of  twelve  strain  gages 
arranged  in  three  bridges  of  four  each.  The  physical  arrangement 
of  the  strain  gages  is  depicted  in  Figure  4,  while  the  electrical 
arrangement  is  presented  in  Figure  5. 

Design  of  the  force  balance  was  complicated  by  the  fact  that  the  bal¬ 
ance  has  to  carry  the  inertia  loads  as  well  as  the  aerodynamic  loads. 
The  moment  inertia  loads  become  quite  large  as  the  oscillation 
frequency  increases.  A  problem  arises  in  trying  to  design  for  the 
sensitivity  required  for  axial  force  measurements  and  still  have  the 
capability  of  transmitting  the  large  moments.  For  oscillating  force 
measurements,  an  additional  problem  of  avoiding  natural  frequencies 
in  the  desired  operating  range  arises.  With  these  problems  in  mind, 
it  was  decided  to  design  the  balance  for  an  oscillating  frequency 
limit  of  15  Hertz  in  anticipation  that  if  the  balance  design  were 
satisfactory  it  could  be  scaled  up  for  later  tests  at  higher  frequen¬ 
cies. 

The  same  carrier  equipment  and  galvanometers  used  for  the  pres¬ 
sure  instrumentation  were  used  to  record  the  force  data.  For  the 
force  data,  only  one  recorder  was  required  since  only  six  channels 
were  needed  for  the  two  balances. 

TESTING  PROCEDURE  AND  ANALYSIS  OF  DATA 


All  testing  was  done  in  the  1.  5 -foot  by  3.  875  foot  wind  tunnel  at  a  dynamic 
pressure  of  28.  205  pounds  per  square  foot  and  an  indicated  airspeed  of 
105  mph.  Actual  test  conditions  (tunnel  temperature  and  pressure)  re¬ 
sulted  in  velocities  ranging  between  106  and  108  mph.  The  Reynolds 
number  based  upon  the  model  chord  length  was  0.  93  x  10®.  Model  oscilla¬ 
ting  frequency  was  varied  from  2  to  15  Hertz  (frequency  parameter  varied 
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from  0.  03  to  0.  3)  in  nine  steps  for  each  test  condition  listed  in  Table  I  . 
Tunnel  wall  corrections  were  not  applied  to  the  data. 


TABLE  I  .  SUMMARY  OF  TEST  CONDITIONS 

Pitch-Axis 

Location,  % 

5 

A<* 

Model 

Chord 

Deg 

Deg 

Pressure 

25 

-0.20 

4.00 

-0.20 

6.  00 

5.  80 

6.00 

13.  83 

6.00 

18.  00 

6.00 

Force 

25 

-0.  35 

6.08 

5.81 

6.08 

13.  56 

6.08  1 

37 

-0.31 

6.08 

5.  84 

6.08 

13.76 

6.08 

50 

-0.  02 

6.08 

6.22 

6.08 

14.25 

6.08 

1 ...  . .  1 

Pressure  Model 


Before  testing  the  pressure  model,  the  pressure  transducers  were  cali¬ 
brated.  This  calibration  was  accomplished  by  connecting  the  reference 
pressure  manifold  to  an  alchohol  manometer  and  applying  various  pressures 
to  the  system.  By  recording  the  manometer  reading  and  taking  an  oscillo¬ 
graphic  record  for  each  pressure,  the  entire  system  was  calibrated  using 
the  alcohol  manometer  as  a  standard.  Since  the  tunnel  dynamic  pressure 
is  determined  by  an  alcohol  manometer  using  the  same  fluid,  the  pressure 
coefficients,  AP/q,  are  also  independent  of  specific  gravity  of  the  alcohol. 
The  reference  pressure  manifold  provided  a  simple  means  of  calibrating 
the  system  and  spot -checking  the  calibration  prior  to  each  run. 

Having  calibrated  the  transducers,  steady  state  data  were  obtained  for 
angles  of  attack  varying  from  -4  degrees  to  30  degrees  in  increments  of 
2  degrees.  Pressure  data  were  recorded  for  each  point.  Data  for  the 
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individual  orifices  were  read  from  the  oscillograph  record  and  punched 
into  IBM  cards  using  an  oscillograph  chart  reader.  The  cards  were  pro¬ 
cessed  by  computer  to  obtain  pressure  coefficients,  normal  force,  and 
pitching  moment.  Normal  forces  were  obtained  by  integrating  the  pressure 
coefficients  using  a  trapezoidal  method.  Both  horizontal  and  vertical 
contributions  of  the  pressures  were  considered  in  the  integration  to  obtain 
moments.  For  purposes  of  computation,  the  trailing  edge  pressure  was 
assumed  to  be  zero. 

Before  conducting  the  oscillating  test,  the  transducers  were  checked  for 
gravitational  effects  on  the  diaphragm  due  to  model  oscillations.  To  accom 
plish  this,  the  pressure  orifices  were  taped  at  the  model  surface.  The 
model  was  oscillated  through  the.  speed  range,  and  the  oscillograph  traces 
were  observed.  With  the  exception  of  the  transducer  located  at  25  percent 
chord  on  the  upper  surface,  there  was  no  indication  of  any  gravitational 
effects.  The  transducer  mentioned  apparently  had  a  loose  connection  or 
part,  for  at  frequencies  above  8  or  10  Hertz,  the  signal  became  erratic. 
Below  this  point  there  was  no  indication  of  any  gravitational  effects. 

In  order  to  check  for  variations  in  tunnel  speed  due  to  model  oscillations, 
a  pressure  transducer  was  connected  to  the  piezometer  ring  just  forward 
of  the  test  section.  The  signal  from  this  transducer  was  observed  with 
the  tunnel  operating,  and  the  model  oscillations  varied  through  the  test 
range.  There  was  no  indication  of  unsteadiness  arising  from  the  oscilla¬ 
ting  model. 

Being  satisfied  that  there  were  no  extraneous  signals  due  to  the  design  of 
the  system  and  operating  conditions,  the  oscillating  tests  listed  in  Table 
I  were  conducted.  These  tests  were  chosen  so  that  the  effects  of  oscilla¬ 
ting  amplitude  and  mean  angle  of  attack  could  be  determined.  The  higher 
mean  angles  were  chosen  so  that  one  would  be  in  the  vicinity  of  the  steady 
state  stall  angle  of  attack  and  the  other  would  be  well  within  the  stall 
region.  For  the  high  mean  angle  of  attack  tests,  the  pressure  transducers 
were  biased  by  applying  a  negative  pressure  to  the  reference  pressure 
manifold.  This  shifted  the  mean  values  so  that  the  oscillating  values 
would  not  exceed  the  deflection  limitations  of  the  recording  equipment. 

In  this  way,  the  sensitivity  of  the  oscillating  signal  was  not  reduced,  as  it 
would  be  if  the  signal  had  been  attenuated  to  keep  the  deflections  down. 

In  order  to  relate  the  pressure  signals  to  the  model  motion,  a  signal  from 
a  potentiometer  attached  to  the  pitch  shaft  was  added  to  both  recorders. 

A  timing  mark  was  applied  to  the  position  signals  to  tie  the  two  records 
together. 

Records  from  the  oscillating  tests  were  read  and  punched  into  cards  in 
the  same  manner  as  the  steady  state  data.  Twenty-four  points  on  each  of 
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three  consecutive  cycles  were  read  for  the  mean  angles  of  attack  for 
which  the  instantaneous  angle  of  attack  did  not  exceed  the  steady  state 
stall  value.  Instantaneous  normal  forces  and  pitching  moments  were 
computed  for  these  points.  These  values  were  then  plotted  against  air¬ 
foil  angular  position.  The  magnitudes  of  the  oscillating  components 
of  normal  force  and  pitching  moment  and  the  phase  relationships  were 
measured  from  these  plots.  For  the  two  higher  mean  angles  of  attack, 
only  one  cycle  was  read  and  computed.  The  nominal  fcpt  value  was  set  on 
the  oscillating  mechanism;  however,  due  to  deflections  in  the  system,  the 
actual  value  of  fax  was  higher  than  the  nominal  value,  and  the  difference 
increased  with  increased  oscillating  frequency.  The  actual  value  of  bflt, 
which  wais  determined  by  measuring  the  travel  of  the  trailing  edge  of  the 
model  at  various  oscillating  speeds,  was  considered  in  the  computation  of 
the  absolute  magnitudes  of  the  normal  force  and  pitching  moment. 

Force  Models 


Calibration  of  the  strain  gage  balances  was  undertaken  after  completion 
of  the  pressure  testing.  The  balances  were  calibrated  in  place  by  substi¬ 
tuting  an  aluminum  plate,  similar  to  the  force  model  spar,  for  the  model. 
The  plate  was  fitted  with  attachments  for  loading  normal  force,  axial 
force,  and  pitching  moments.  The  balances  were  loaded  in  increments  of 
the  primary  loads  to  values  slightly  higher  than  the  anticipated  test  loads. 
In  addition  to  the  primary  loada  combined  loads  were  applied  to  determine 
balance  interactions.  No  second-order  interactions  of  any  consequence 
were  detected  from  the  results  of  the  combined  loads  tests.  Interaction 
equations  were  written  for  each  balance  to  account  for  the  primary  inter¬ 
actions.  Sensitivity  of  the  normal  force  and  pitching  moment  was  very 
good,  but  axial  force  was  less  sensitive  than  was  desirable. 

The  force  model  ifrith  the  pitch  axis  at  25  percent  chord  was  installed  to 
check  out  the  system,  oscillating  mechanism,  balances,  and  instrumenta¬ 
tion.  At  low  oscillating  frequencies,  the  oscillograph  traces  were  quite 
smooth.  As  the  frequency  of  oscillation  was  increased,  a  higher  fre¬ 
quency  disturbance  appeared  on  the  traces.  The  magnitude  of  this  distur¬ 
bance  increased  as  oscillating  frequency  increased  and  reached  alarming 
proportions  at  the  highest  frequencies.  This  disturbance  was  found  to  be 
arising  from  the  excitation  of  the  natural  frequency  of  the  model-balance 
system.  The  natural  frequency  in  the  pitch  mode  was  approximately  130 
Hertz  while  that  of  the  axial  force  mode  was  about  30  Hertz.  Considerable 
effort  was  expended  in  trying  to  reduce  the  magnitude  of  this  disturbance 
by  damping  and  by  increasing  the  natural  frequency  so  that  it  would  not  be 
excited  as  readily.  Some  damping  was  achieved  by  using  heavier  supports 
and  a  bearing  on  the  pitch  shaft  as  close  to  the  balance  as  possible.  The 
only  way  the  natural  frequency  could  be  increased  was  to  decrease  the 
model  mass  or  increase  the  stiffness  of  the  balances.  The  weight  of  the 
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model  (2.  63  pounds)  could  not  be  decreased  appreciably.  Since  the  sensi¬ 
tivity  of  the  axial  force  balance  was  already  less  than  that  desired,  the 
balance  could  not  be  stiffened.  Therefore,  it  was  decided  to  proceed  with 
the  tests  and  depend  upon  harmonic  analysis  of  the  data  to  provide  meaning¬ 
ful  results.  Fortunately  when  the  model  was  oscillated  with  the  wind  tunnel 
running,  there  was  an  appreciable  damping  of  the  natural  frequency  signal. 

Since  the  force  balances  measure  the  inertia  loads  as  wel1  as  the  aero¬ 
dynamic  forces,  it  was  necessary  to  determine  these  loads  and  subtract 
them  from  the  total  values  to  obtain  the  desired  aerodynamic  character¬ 
istics.  An  attempt  was  made  to  oscillate  the  model  in  a  vacuum  to 
measure  the  moments  of  inertia  independent  of  the  virtual  mass  effect  of 
the  air.  The  model  was  oscillated  in  a  tank  at  atmospheric  pressure  and 
then  with  the  tank  evacuated  to  a  pressure  of  27  inches  of  mercury  below 
atmospheric.  Leakage  around  the  shaft  prevented  the  attainment  of  a 
higher  vacuum.  There  was  no  measurable  difference  between  the  results 
of  these  two  tests.  Therefore,  the  virtual  mass  effect  of  the  air  was 
assumed  to  be  negligible,  and  the  inertia  loads  were  determined  by  oscilla¬ 
ting  the  model  in  still  air. 

The  procedure  of  testing  the  three  force  models  was  the  same  for  each  one. 
First  the  models  were  mounted  in  the  tunnel  on  the  oscillating  mechanism. 
Known  loads  were  applied  to  the  model  to  check  the  balance  calibrations. 
Then  steady  state  wind-on  runs  were  made  varying  the  angle  of  attack  from 
-4  degrees  to  30  degrees,  data  being  recorded  for  2 -degree  increments 
through  the  stall  angle  and  4 -degree  increments  above  stall.  After  obtain¬ 
ing  the  steady  state  data,  the  desired  mean  angle  was  set  on  the  oscillator 
and  the  model  was  oscillated  through  the  frequency  range  with  the  wind  off 
to  obtain  the  inertia  loads.  Data  were  recorded  at  nine  different  oscilla¬ 
ting  frequencies.  Immediately  after  recording  the  inertia  loads,  the 
tunnel  was  brought  up  to  speed,  the  model  was  oscillated  at  the  same  fre¬ 
quencies  employed  for  the  inertia  loads,  and  data  were  recorded. 

Oscillograph  records  were  read  and  punched  into  cards  in  the  same  manner 
as  the  pressure  data.  In  this  case,  only  seven  channels  were  required, 
three  for  each  balance  and  one  for  a  position  trace.  The  data  were  then 
processed  by  computer.  Forces  and  moments  were  computed  for  the 
balances  using  the  interaction  equations.  These  results  for  the  two  balances 
were  then  averaged  to  obtain  the  model  normal  force,  axial  force,  and 
pitching  moment  (in  all  cases  pitching  moment  is  aboul  the  pitch  axis) 
acting  on  the  model.  A  24-point  harmonic  analysis  was  then  performed  on 
both  inertia  data  and  wind-on  data.  The  results  of  the  inertia  analysis 
were  subracted  from  the  wind -on  results  to  obtain  the  final  aerodynamic 
coefficients  and  phase  relationships. 

As  was  mentioned  in  the  discussion  on  the  pressure  model,  actual  tip 
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values  were  determined  by  measurements.  The  deflection  of  the  model 
was  also  measured  with  a  pitching  moment  applied.  The  correction  to 
the  value  obtained  from  the  two  methods  exhibited  excellent  agreement. 
All  of  the  nominal  &at  val-es  were  corrected  for  deflections  which  were  a 
function  of  the  oscillating  frequency. 
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RESULTS  AND  COMPARISON  OF  DATA 


Result*  of  the  oscillating  tests  conducted  at  low  mean  angles  of  attack 
(approximately  0  and  6  degrees)  are  presented  in  coefficient  form  defined 
as  follows: 


dC 


N 


da 


Ini 


1/2  pV  S  &a 


and 

dCM  ,  IMI 

da  l/2pV2Sc  & 

where  £a  is  measured  in  radians. 

For  mean  angles  where  instantaneous  angles  of  attack  exceed  the  steady 
state  stall  angle  of  attack,  representative  data  are  presented  as  instan¬ 
taneous  coefficients  versus  instantaneous  angle  of  attack. 

STEADY  STATE  RESULTS 


Steady  state  Cn  and  Cj^  variations  with  a  are  presented  in  Figure  6.  The 
slopes  of  the  Cn  versus  a  curves  for  all  four  models  show  excellent  agree¬ 
ment.  However,  there  is  considerable  variation  of  the  curves  in  the  region 
near  stall.  This  variation  points  up  some  of  the  problems  associated  with 
airfoil  stall.  All  of  these  models  were  cast  from  the  same  mold  and 
tested  in  the  same  facility.  Strictly  speaking,  the  steady  state  Cn  versus 
a  curves  should  be  the  same.  It  is  assumed  that  slight  irregularities  at 
the  leading  edge  of  the  surface  cause  the  stall  to  be  precipitated  differ¬ 
ently  on  each  of  the  four  models.  In  addition  to  the  problem  of  initial 
stall  is  the  problem  of  the  unsteady  and  irregular  nature  of  the  flow  over 
the  model  after  it  has  stalled.  Oscillograph  records  taken  in  the  stall 
region  show  rapid  fluctuations  of  the  order  of  30-50  percent  of  the  maxi¬ 
mum  recorded  normal  force.  Data  presented  in  Figure  6  for  the  stall 
region  are  average  values  of  these  fluctuations.  The  discrepancies  noted 
in  this  relatively  simple  case  of  steady  state  stall  are  emphasized  to 
illustrate  the  problem  of  determining  maximum  Cn*  stall  angle  of  attack, 
and  aerodynamic  characteristics  after  the  inception  of  stall. 
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PRESSURE  MODEL  RESULTS 


Normal  Force  and  Pitching  Moment  Coefficients 

The  magnitude  and  phase  of  the  oscillatory  normal  force  and  pitching 
moment  coefficients  for  low  mean  angles  of  attack  are  presented  in  Figures 
7  and  8  as  functions  of  reduced  frequency.  For  purposes  of  comparison, 
the  theoretical  normal  forces  and  pitching  moments  calculated  from 
Theodorsen's  equations^  are  shown  along  with  the  measured  data.  Meas¬ 
ured  data  are  presented  for  two  values  of  mean  angle  of  attack  and  two 
values  of  oscillating  amplitude.  All  of  the  data  shows  excellent  agreement 
with  theory  except  for  the  moment  phase  angles.  The  oscillating  ampli¬ 
tude  of  the  pitching  moment  is  quite  small  with  the  pitch-axis  located  at 
25  percent  chord/  especially  at  low  oscillating  frequencies.  This  presents 
a  problem  in  trying  to  measure  phase  angles  with  any  degree  of  accuracy 
and  probably  accounts  in  large  measure  for  the  scatter  in  pitching  moment 
phase  data  and  the  discrepancy  with  theory.  Oscillatory  amplitudes  of  nor¬ 
mal  force  and  pitching  moment  coefficients  and  phase  relationships  appear 
to  be  independent  of  mean  angle  of  attack  and  amplitude  of  oscillation  as 
long  as  the  instantaneous  value  of  angle  of  attack  does  not  exceed  the 
steady  state  stall  value. 

Instantaneous  normal  force  and  pitching  moment  coefficients  are  presented 
in  Figures  9  and  10  as  functions  of  instantaneous  angle  of  attack  for  two 
values  of  mean  angle  of  attack  at  representative  test  values  of  reduced 
frequency.  Figure  9  presents  data  for  mean  angle  of  attack  a  of  13.  80 
degrees,  which  is  close  to  the  value  of  a  for  maximum  formal  force. 

Figure  10  is  for  a  value  of  a  of  18  degrees,  which  is  well  above  the  steady 
state  stall  angle  of  attack.  Steady  state  normal  force  and  pitching  moment 
are  included  in  Figures  9  and  10  for  comparison  with  the  oscillatory  values. 
In  Figure  9  it  is  observed  that  at  the  lowest  values  of  reduced  frequency, 
the  model  stalls  at  an  angle  of  attack  slightly  greater  than  the  steady  state 
stall  angle.  The  pitching  moment  increases  in  magnitude  to  a  large 
negative  value  that  is  considerably  greater  than  the  steady  state  stall 
value.  With  increasing  values  of  k  and  reduced  frequency,  the  angle  of 
stall  for  the  normal  force  is  delayed  until  at  the  highest  value  of  k,  there 
is  essentially  no  indication  of  stall.  A  maximum  value  of  Cn  of  approxi¬ 
mately  1.  8  is  obtained  at  the  highest  value  of  k.  This  is  an  increase  of 
about  40  percent  over  the  steady  state  value  of  maximum  CN.  The  increase 
in  negative  at  stall  decreases  as  k  increases  and  approaches  the 
unstalled  condition  at  the  highest  frequency.  Figure  10  indicates  much 
the  same  trends  as  Figure  9.  As  k  increases,  the  angle  for  CN  stall 
increases.  At  the  highest  value  of  k,  the  maximum  value  of  Cn  is  not 
reached  until  after  the  maximum  value  of  a  has  been  obtained.  The  maxi¬ 
mum  value  of  CN  (2.  29)  for  the  highest  oscillating  frequency  exceeds  the 
steady  state  value  by  approximately  80  percent.  The  pitching  moment 


15 


stalls  somewhat  sooner  than  the  normal  force  in  all  cases.  As  k  increases, 
the  pitching  moment  at  stall  increases  in  magnitude  until  a  value  of  -0.  33 
is  reached  at  k  =  0.  191.  With  further  increase  in  k,  this  large  negative 
moment  appears  to  decrease  in  a  manner  similar  to  the  case  with  S  - 
13.  80  degrees. 


Instantaneous  Pressure  Distributions 


Instantaneous  pressure  distributions  are  presented  in  Figures  11  through 
16  for  some  representative  conditions.  Figures  11  and  12  present  pres¬ 
sure  distributions  for  a  =  5.  80  degrees  at  the  extremes  of  the  frequency 
range.  Data  are  presented  for  an  instantaneous  angle  of  attack  to  compare 
the  instantaneous  pressure  distributions  for  that  portion  of  the  cycle  where 
a  is  increasing  with  the  pressure  distribution  for  the  same  angle  when  a 
is  decreasing.  Figure  11  is  representative  of  the  low-frequency  end  of  the 
data  summarized  in  Figures  7  and  8.  From  Figure  11,  it  is  seen  that 
there  is  very  little  difference  between  increasing  a  and  decreasing  a. 

This  is  in  agreement  with  Figures  7  and  8,  which  indicate  a  very  small 
normal  force  phase  angle  at  low  frequency  and  very  small  pitching  mo¬ 
ments.  For  the  high-frequency  data  of  Figure  12,  there  is  a  noticeable 
difference  in  the  normal  force  (area  under  curve)  and  pitching  moment 
(area  distribution)  between  increasing  and  decreasing  a.  The  normal  force 
is  greater  for  a  increasing  than  for  a  decreasing.  This  is  indicative  of  the 
leading  phase  angle  shown  in  Figure  7.  The  pitching  moment,  essentially 
zero  for  a  increasing,  becomes  positive  with  decreasing  a,  which  is  indi¬ 
cative  of  the  pitching  moment  phase  lag  shown  in  Figure  8.  The  pressure 
distributions  also  indicate  a  decrease  in  Cjj  magnitude  between  the  lowest 
and  highest  reduced  frequencies  (Figure  7  and  8  respectively),  but  this  is 
not  so  obvious  without  overlaying  one  figure  with  the  other. 

Figures  13  through  15  present  instantaneous  pressure  distributions  for 
selected  values  of  a  for  low,  medium,  and  high  test  frequencies  respec¬ 
tively.  The  mean  angle  a  for  these  figures  was  13.  80  degrees.  Data  are 
presented  for  instantaneous  values  of  a  near  the  mean  and  near  the  maxi¬ 
mum  values.  Figure  13  indicates  that  the  model  stalls  before  the  maxi¬ 
mum  angle  of  attack  is  reached  at  the  low  oscillating  frequency.  For  the 
midfrequency  value  presented  in  Figure  14,  the  model  stalls  at  or  very 
near  the  maximum  angle  of  attack,  as  indicated  by  the  stalled  condition 
for  decreasing  a.  Figure  15  indicates  that  the  model  is  essentially 
unstalled.  These  figures  aid  in  interpreting  the  results  in  Figure  9. 

Figure  16  presents  some  very  unusual  pressure  distributions  obtained  from 
the  highest  oscillating  frequency  test  for  a  mean  angle  of  attack  of  18.  00 
degrees.  Data  are  presented  for  several  instantaneous  angles  of  attack 
over  the  positive  half  cycle  of  &a.  These  pressure  distributions  are 
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directly  related  to  the  instantaneous  force  and  moment  data  presented  in 
Figure  10  for  the  highest  frequency.  It  is  interesting  to  note  in  Figure  16 
that  Cjj  increases  up  to  and  beyond  maximum  angle  of  attack  even  though 
ihe  peak  suction  pressures  near  the  leading  edge  drop  off.  The  relatively 
large  negative  pressure  coefficients  existing  over  the  airfoil  upper  sur¬ 
face  aft  of  the  25  percent  chord  are  responsible  for  the  large  negative 
moments  shown  in  Figure  10.  Also  of  interest  is  the  flat  pressure  distri¬ 
bution  over  the  first  20  percent  of  the  lower  surface  of  the  airfoil. 

Figure  17  presents  instantaneous  pressure  distributions  for  the  lower  sur¬ 
face  leading  edge  of  the  airfoil  for  various  instantaneous  angles  of  attack, 
with  increasing  and  decreasing  a  for  oscillations  about  a  mean  angle  of 
5.80  degrees.  The  irregular  nature  of  the  distributions  presented  here 
was  not  noted  for  any  of  the  other  test  conditions.  These  irregularities 
may  be  indicative  of  a  vortex  formation  or  some  other  flow  peculiarity. 

In  the  future,  it  may  be  desirable  to  do  some  sort  of  flow  visualization  on 
the  airfoil  oscillating  at  this  condition. 

Representative  oressure  distributions  were  chosen  for  presentation  in 
this  report.  However,  since  there  is  little  instantaneous  experimental 
pressure  distribution  data  available,  it  was  felt  that  tabulated  coefficients 
should  be  presented  for  one  cycle  of  each  test  condition.  These  coeffi¬ 
cients  are  listed  in  Appendix  111. 

FORCE  TEST  RESULTS 


Normal  Force  and  Pitching  Moment  Coefficients 

Magnitude  and  phase  of  the  oscillatory  normal  force  and  pitching  moment 
coefficients  for  the  model  oscillating  about  the  25  percent  chord  axis  at 
low  mean  angles  of  attack  are  presented  in  Figure  18  and  19.  As  was  the 
case  with  the  pressure  model,  the  experimental  data  agree  well  with 
theory  except  for  the  moment  phase  angles.  The  moment  phase  data  show 
reasonable  agreement  at  the  highest  values  of  k,  but  diverge  from  theory 
at  the  lower  k  values, tending  toward  zero  phase  angle  at  k  equal  zero. 
Admittedly,  the  accuracy  of  the  pitching  moment  data  in  this  region 
leaves  much  to  be  desired,  but  there  appears  to  be  a  definite  trend  in 
both  sets  of  data.  An  error  of  just  two  counts  of  pitching  moment  (.  002) 
can  produce  an  error  of  greater  than  10  degrees  phase  angle  in  the  low  k 
region  due  to  the  small  magnitudes  of  the  moment.  But  it  would  be  ex¬ 
pected  that  errors  would  produce  scatter  and  not  such  a  noticeable  trend. 

Oscillatory  coefficients  for  the  model  oscillating  about  ihe  37  percent 
chord  are  presented  in  Figures  20  and  21.  There  is  excellent  agreement 
between  measured  values  and  theory  except  for  the  moment  phase  angles 
at  the  higher  k  values.  Here  there  is  a  tendency  to  diverge  from  theory, 
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with  measured  phase  angles  being  less  than  the  theoretical  values. 

Data  presented  in  Figures  22  and  23  for  oscillations  about  the  50  percent 
chord  also  exhibit  excellent  agreement  with  theory.  The  measured 
moment  values  agree  very  well  with  theory,  but  there  still  appears  to  be  a 
tendency  to  diverge  at  the  higher  reduced  frequencies,  as  was  evidenced 
in  the  case  of  the  37  percent  pitch  axis. 


Drag  Coefficients 

As  was  mentioned  previously,  the  d£ag  balance  was  not  as  sensitive  as  was 
desired.  An  error  in  reading  the  oscillograph  record  of  .  01  in.  was  equiva¬ 
lent  to  .0030  (30  counts)  in  drag  coefficient.  For  the  0012  airfoil  tested, 
this  is  approximately  50  percent  of  the  minimum  drag.  This  low  sensi¬ 
tivity  made  it  difficult  to  obtain  reasonable  steady  state  data  at  low  angles 
of  attack  where  the  drag  is  quite  low.  In  addition  to  the  sensitivity  problem, 
a  temperature  problem  existed  on  the  drag  balances.  This  apparently  comes 
from  the  wide  spacing  of  the  strain  gages  in  the  drag  bridge  (see  Figure  4). 
The  carrier  voltage  heats  the  gages.  The  gages  are  then  cooled  by  air  cir¬ 
culating  around  them  when  the  tunnel  is  operating.  Some  air  can  circulate 
since  there  is  a  small  gap  between  the  model  and  the  tunnel  wall.  This 
circulation  of  air  causes  unequal  cooling  at  the  drag  strain  gages.  This 
cooling  is  negligible  on  the  pitching  moment  and  normal  force  bridges. 

This  temperature  shift  varied  on  the  three  models  tested,  being  greatest 
on  the  25  percent  chord  pitch  axis  model  and  negligible  on  the  50  percent 
chord  pitch  axis  model.  A  test  technique  eliminated  most  of  this  tempera¬ 
ture  drift.  The  tunnel  was  brought  up  to  speed  and  allowed  to  run  until  the 
drag  values  stabilized.  Then  the  tunnel  was  shut  down,  wind-off  zeros  were 
taken,  and  the  tunnel  was  immediately  started  again.  The  small  remaining 
error  (arising  from  the  finite  time  required  to  stop  all  flow  in  the  tunnel) 
was  sufficient  to  cause  problems  due  to  the  low  sensitivity  of  the  balance. 
This  problem  was  especially  aggravating  in  trying  to  obtain  steady  state 
and  mean  values  of  drag.  It  did  not  affect  the  oscillating  values,  since 
there  was  negligible  drift  during  the  time  required  to  record  the  oscilla¬ 
ting  data.  However,  obtaining  the  oscillating  data  is  complicated  by  the 
excitation  of  the  natural  frequencies  of  the  balances. 

In  spite  of  the  difficulties,  some  drag  results  were  obtained  and  are  pre¬ 
sented  in  Figures  24  and  25.  Drag  coefficients  were  obtained  from  the 
axial  force  data  by  means  of  the  following  relationship: 

C_>  =  C .  cos  a  +  C.T  sin  a 
DA  N 

Unlike  the  normal  force  prior  to  stall,  drag  variation  is  not  linear  with 
angle  of  attack.  Therefore,  axial  force  is  not  linear  with  angle  of  attack. 

The  axial  force  data  obtained  from  the  harmonic  analysis  indicate  that  the 
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first  three  harmonics  are  significant.  Higher  harmonics  are  negligible, 
except  those  arising  from  the  natural  frequencies  which  are  unwanted.  The 
first  and  second  harmonics  contain  most  of  the  data  of  interest.  The  third 
harmonic  was  small  in  comparison  to  the  first  two  and  was  essentially  con¬ 
stant.  Since  it  was  constant  with  k  and  relatively  small,  neglecting  it  does 
not  alter  the  drag  picture.  Data  presented  in  Figure  24  are  for  a  mean 
angle  of  attack  of  zero  degrees.  Figure  25  presents  the  dat8  for  a  mean 
angle  of  approximately  6  degrees.  Oscillating  amplitude  is  a  nominal  6 
degrees  in  each  case.  Oscillating  drag  is  entirely  different  for  the  two 
mean  angle  conditions.  In  the  case  of  zero  mean  angle,  drag  increases 
with  both  positive  and  negative  angles  of  attack.  Since  the  drag  variation 
with  angle  of  attack  is  fairly  flat  in  this  low  a  region,  the  drag  results 
would  be  expected  to  be  essentially  a  double -frequency  sinusoid.  For  a 
nominal  mean  angle  of  6  degrees,  the  drag  decreases  to  a  minimum  value 
for  -  A®  and  increases  to  a  maximum  for  +  A®.  The  drag  curve  over  this 
range  is  nonlinear,  with  CD  increasing  more  rapidly  with  a  at  the  higher 
angles  of  attack.  Therefore,  it  would  be  expected  that  drag  would  be  a 
nonsinusoidal  curve  of  the  oscillating  frequency. 

Results  of  the  harmonic  analysis  of  the  zero  mean  angle  of  attack  data  show 
very  small  first  harmonic  and  third  harmonic  contributions  to  the  drag. 

The  magnitude  of  the  first  harmonic  data  is  probably  due  to  the  fact  that 
the  mean  angle  of  attack  was  not  exactly  zero  degrees.  The  second  har  - 
monic  data  are  presented  along  with  the  mean  drag  coefficient  in  Figure  24 
as  a  function  of  oscillatory  frequency  parameter.  The  mean  value  of  drag 
is  seen  to  increase  rapidly  with  k.  The  second  harmonic  is  constant  up  to 
a  value  of  k  =  0.  2,  then  increases  sharply.  The  sharp  increase  in  the 
second  harmonic  is  probably  due  to  the  excitation  of  the  balance  natural 
frequency,  since  the  natural  frequency  in  the  drag  direction  is  approximately 
30  Hertz,  which  is  the  second  harmonic  of  the  oscillating  frequency  of 
15  Hertz  at  k  =  0.  3. 

Harmonic  analysis  of  the  nominal  6-degree  mean  angle  of  attack  data  re¬ 
sults  in  large  values  of  the  first  harmonic,  sizeable  values  of  the  second 
harmonic,  and  small  values  of  the  third  harmonic.  The  second  and  thiro 
harmonics  probably  arise  from  the  nonlinearity  of  the  drag  curve.  Figure 
25  presents  the  mean  drag  coefficient  and  the  first  two  harmonics  as 
functions  of  the  frequency  parameter.  As  was  the  case  for  a  =  0,  the  mean 
value  of  drag  increases  rapidly  with  increasing  k.  The  amplitudes  of  the 
first  and  second  harmonics  decrease  with  increasing  k. 

Data  for  both  mean  angles  of  attack  exhibit  the  rapid  increase  of  C  with 
increasing  frequency  of  oscillation.  For  low  values  of  k,  one  woura  expect 
the  oscillating  drag  variation  to  follow  the  steady  state  variation  with 
angle  of  attack.  Mean  drag  coefficients  at  these  low  oscillation  frequencies 
are  lower  than  anticipated.  Some  of  this  discrepancy  may  result  from  the 
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drag  balance  temperature  problem  discussed  previously,  but  it  is  not  con¬ 
ceivable  that  all  of  the  discrepancy  could  be  arising  from  this  source. 

For  the  steady  state  condition,  the  drag  is  a  minimum  at  or  =  0  degrees  and 
the  area  of  the  model  surface  experiencing  turbulent  flow  is  a  minimum. 

As  a  increases,  the  area  of  turbulent  flow  increases.  With  the  model 
oscillating,  viscous  effects  may  keep  the  turbulent  area  from  returning  to 
the  minimum  condition,  thereby  keeping  the  drag  from  obtaining  the  mini¬ 
mum  value.  As  k  increases,  the  low  drag  diverges  from  the  steady  state 
minimum.  This  would  account  for  an  increase  in  the  mean  value  of  drag. 

It  would  also  account  for  the  more  rapid  rise  4n  for  5  =  6  degrees, 
since  the  steady  state  drag  increase  from  a  =  6  degrees  to  a  =  12  degrees 
is  roughly  four  times  as  great  as  the  drag  rise  from  0  to  6  degrees.  If 
this  is  really  the  case,  then  the  oscillating  component  of  drag  should  de¬ 
crease  with  increasing  k.  This  is  borne  out  by  the  decrease  in  the  first 
harmonic  amplitude  for  the  model  oscillating  about  a  =  6  degrees.  How¬ 
ever,  there  is  no  evidence  of  a  decrease  in  the  amplitude  of  the  second 
harmonic  for  the  or  =  0  condition.  This  may  be  because  the  change  in  ampli¬ 
tude  for  this  condition  would  be  so  small  that  it  is  masked  by  the  scatter 
in  the  data. 

In  the  data  presented  for  a  3  0  degrees,  there  is  no  indication  of  any  changes 
in  the  drag  coefficients  due  to  change  in  pitch-axis  location.  For  5  =  6 
degrees,  there  are  considerable  changes  with  pitch  axis,  especially  in  the 
first  harmonic  values.  The  mean  value  of  drag  tends  to  increase  as  the 
pitch  axis  is  moved  rearward  and  the  oscillatory  amplitude  (first  harmonic) 
decreases.  The  large  change  in  the  first  harmonic  may  be  due  to  an  in¬ 
creased  camber  effect  as  the  pitch  axis  is  moved  aft.  The  leading  edge  of 
the  model  has  increased  travel  as  the  pitch  axis  is  moved  rearward,  caus¬ 
ing  an  increase  in  induced  camber.  The  effect  of  the  camber  is  to  shift 
the  drag  curve  so  that  the  minimum  drag  occurs  at  a  higher  angle  of  attack. 
As  the  drag  curve  is  shifted,  the  oscillation  takes  place  over  a  flatter  re¬ 
gion  of  the  curve,  thereby  reducing  the  difference  between  minimum  and 
maximum  values. 

The  second  harmonic  data  presented  in  Figure  25  for  a  =  6  degrees  show 
a  decrease  in  amplitude  with  increasing  k  (varying  somewhat  with  pitch- 
axis  location).  If  the  first  harmonic  (amplitude)  is  decreasing  with  k,  then 
a  decrease  in  the  second  harmonic  may  result  from  an  accompanying  in¬ 
crease  in  linearity.  It  should  be  noted  that  speaking  of  amplitudes  of  the 
harmonics  is  not  the  same  as  amplitude  of  oscillating  drag,  since  the  drag 
is  composed  of  the  sum  of  the  harmonics. 

In  spite  of  the  problems  associated  with  obtaining  the  drag  data  and  the 
scatter  of  the  results,  it  is  felt  that  the  drag  trends  are  quite  pronounced; 
the  consistency  obtained  from  the  three  models  verifies  these  trends.  It 
is  felt  that  the  large  mean  drag  increases  reported  here  would  warrant 
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further  investigation.  A  drag  balance  of  the  type  used  may  be  feasible  if 
the  sensitivity  can  be  increased.  It  may  be  possible  to  use  semiconductor 
strain  gages  on  the  balance  which  give  an  order  of  magnitude  increase  in 
sensitivity  over  the  foil  gages  presently  employed.  If  this  were  possible, 
then  the  beams  may  be  thickened,  pushing  the  natural  frequency  up,  and 
the  sensitivity  may  still  be  increased  by  a  factor  of  two  or  three. 

Instantaneous  Normal  Force  and  Moment  Coefficients 

Instantaneous  normal  force  and  moment  coefficients  are  presented  in 
Figures  26  and  27  as  functions  of  instantaneous  angle  of  attack  for  the  30 
percent  chord  model  oscillating  about  a  mean  angle  of  attack  of  6.  22 
degrees.  Figures  26  and  27  are  for  low  oscillating  frequency  and  high 
oscillating  frequency  respectively.  These  figures,  along  with  Figures  18 
through  23,  summarize  the  normal  force  and  moment  ch  *  acteristics  for 
the  model  oscillating  in  the  linear  portion  of  the  angle  of  attack  range.  The 
50  percent  chord  axis  was  chosen  for  these  figures  because  the  greater 
slope  of  the  pitching  moment  versus  angle  of  attack  curve  aids  in  the  pic¬ 
torial  presentation.  Data  for  the  other  models  could  have  been  presented 
and  the  conclusions  would  not  be  altered.  The  experimental  data  in  Figures 
26  and  27  are  the  fundamental  harmonics  obtained  from  the  results  of  the 
harmonic  analysis.  Results  of  the  harmonic  analysis  indicate  that  there 
was  very  little  deviation  from  pure  sinusoidal  motion.  The  average  amount 
of  second  harmonic  present  in  pitching  moment,  normal  force,  and  motion 
was  equal  to  or  less  than  1  percent  of  the  fundamental  amplitude.  This 
justifies  the  use  of  the  fundamental  harmonic  for  data  presentation  and 
accounts  for  the  smooth  curves  presented.  Theoretical  values  are  pre¬ 
sented  in  the  figures  for  comparison  with  experimental  results. 

Agreement  betv/een  theory  and  experimental  results  is  excellent  as  far  as 
shape,  magnitude,  and  direction  of  traverse  are  concerned;  but  at  the  higher 
frequency,  there  is  a  noticeable  displacement  between  the  curves.  The 
direction  of  traverse  and  the  oblateness  of  the  curve  are  functions  of  the 
phase  relationship.  The  direction  of  traverse  changes  or  the  normal 
force  for  the  low- and  high-frequency  curves.  This  agrees  with  the  data 
of  Figure  22  which  show  the  phase  angle  to  be  about  equal  for  the  two 
frequencies  but  of  different  sign.  The  decrease  in  amplitude  for  the  two 
frequencies  as  indicated  in  Figure  22  shows  up  as  a  tilting  of  the  curve 
in  Figure  27.  The  agreement  between  theory  and  experiment  regarding 
size,  shape,  and  direction  of  traverse  relates  to  the  oscillatory  components. 
However,  the  displacement  of  the  curves  indicates  a  discrepancy  between 
the  mean  values  predicted  by  theory  and  the  results  of  the  tests.  The  mean 
values  of  normal  force  and  moment  predicted  by  theory ^  are: 
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CN  =  2  it  a  C(k)* 

CM  =  0/2  +  a)  it  a  C(k) 

where  a  is  a  constant  for  a  given  model 

C(k)  is  Theodorsen's  complex  circulation  function,  F(k)  +  i  G(k) 

For  k  =  0,  F(k)  is  equal  to  1,  G(k)  equals  zero,  and  steady  state  theory  re¬ 
sults.  At  k  >  0,  C(k)  is  always  less  than  1,  resulting  in  a  decrease  of 
mean  values  with  increasing  oscillating  frequency.  Test  results  give  no 
evidence  of  this  decrease  of  mean  values.  Some  of  the  test  results  for 
normal  force  mean  values  are  presented  along  with  theory  in  Figure  28. 

Instantaneous  normal  force  and  pitching  moment  coefficients  are  presented 
in  Figures  29  through  31  for  the  25,  37,  and  50  percent  pitch-axis  models, 
respectively,  as  functions  of  instantaneous  angle  of  attack.  The  mean 
angle  for  each  case  is  in  the  proximity  of  the  steady  state  stall  angle  of 
attack.  The  normal  force  curves  are  quite  similar  for  the  three  models 
and  show  the  same  trends  as  the  pressure  model  data.  That  is,  normal 
force  stall  is  delayed  as  the  oscillating  frequency  is  increased  until  at 
high  frequency,  there  is  essentially  no  indication  of  stall.  While  the 
curves  and  the  trends  are  quite  similar,  there  are  some  differences  in  the 
individual  shapes.  These  differences  probably  arise  from  the  differences 
in  the  steady  state  stall  characteristics  of  the  models  and  the  differences 
between  steady  state  stall  angle  of  attack  and  oscillating  mean  angle.  In 
other  words,  if  the  same  model  were  oscillated  at  slightly  different  mean 
angles  close  to  the  steady  state  stall  angle,  the  resulting  instantaneous 
curves  would  have  slightly  different  shapes.  However,  the  general  trend 
with  k  should  remain  the  same.  From  the  oscillograph  records,  it  was 
observed  that  there  was  even  a  slight  variation  in  instantaneous  forces 
and  moments  from  cycle  to  cycle.  The  data  presented  in  the  figures  of  this 
report  are  for  a  representative  cycle. 

Instantaneous  pitching  moments  exhibit  considerable  variation  with  pitch- 
axis  location.  This  would  be  expected  due  to  the  large  differences  in  the 
steady  state  pitching  moment  for  the  different  pitch  axis.  However,  there 
are  also  differences  in  the  trends  with  changes  in  frequency  parameter. 

For  the  pitch  axis  at  the  25  percent  chord  location,  the  pitching  moment 
experiences  a  sharp  negative  increase  in  magnitude  to  a  value  of  -0.  36  at 
k  =  .  103.  As  k  increases,  this  magnitude  is  reduced.  The  pressure  model 


C  and  C  are  used  interchangeably,  since  they  are  very  nearly  equal  at 
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low  angles  of  attack. 
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exhibited  a  similar  increase  in  negative  magnitude  to  -0.  33,  which  is 
almost  an  identical  magnitude,  under  slightly  different  test  conditions. 

The  exact  frequency  at  which  this  increase  in  magnitude  occurs  is  probably 
dependent  upon  the  nature  of  the  steady  state  stall  and  the  mean  angle  at 
which  the  model  is  oscillating.  An  increase  of  reduced  magnitude  r  o 
occurs  for  the  37  percent  pitch-axis  model  at  the  same  k  value.  1.' 
again,  negative  peak  in  Cm  diminishes  as  k  increases.  There  is  no  indica¬ 
tion  of  this  sharp  change  when  the  pitch  axis  is  located  at  the  50  percent 
chord  point.  Here  the  minimum  pitching  moment  is  about  the  same  for  all 
values  of  k.  Thus,  the  sharp  change  in  pitching  moment  appears  to  be  a 
function  of  reduced  frequency  and  pitch -axis  location. 

COMPARISON  OF  DATA 
Comparison  With  Theory 

Theoretical  results  based  on  the  thin  airfoil  incompressible  theory  of 
Theodorsen^  have  been  presented  along  with  the  experimental  results  where 
applicable.  The  agreement  between  theory  and  experiment  is  excellent  for 
normal  force  coefficients  (both  phase  and  amplitude).  Pitching  moment 
amplitude  agrees  well  with  theory,  but  there  are  some  discrepancies  in 
phase  angles.  The  phase  angle  for  the  25  percent  pitch-axis  force  model 
tends  to  approach  zero  degrees  as  k  goes  to  zero  rather  than  the  value  of 
270  degrees  predicted  theoretically.  This  trend  has  also  been  noted  by 
Wyss  and  Herrera*  for  different  airfoil  profiles  oscillating  about  the  quar¬ 
ter  chord.  There  is  no  evidence  of  this  trend  in  the  pressure  data  results, 
but  then  it  was  not  possible  to  determine  the  phase  angles  at  k  values  be¬ 
low  0. 1  with  any  degree  of  accuracy  from  the  test  results.  Theoretically, 
the  moment  amplitude  goes  to  zero  as  k  approaches  zero  (steady  state 
theory).  However,  the  experimental  results  indicate  finite  values  of  C,, 
for  a  other  than  zero.  Therefore,  one  would  expect  an  oscillatory  ampli¬ 
tude  greater  than  zero  at  very  low  k  values,  and  since  the  moment  is  finite, 
as  the  a  range  is  traversed  very  slowly,  the  phase  angle  should  be  essen¬ 
tially  zero  degrees.  Phase  angle  data  for  the  37  percent  pitch -axis  model 
diverge  from  the  theoretical  values  at  higher  k  values.  No  reasonable 
explanation  for  this  divergence  is  now  available,  since  the  data  for  both 
the  25  percent  axis  model  and  the  50  percent  axis  model  show  good  agree¬ 
ment  with  theory  at  the  higher  values  of  k. 

Comparison  of  Pressure  Model  and  Force  Model  Results 

Oscillatory  normal  force  and  pitching  moment  data  obtained  by  pressure 
measurements  and  those  obtained  by  direct  force  measurements  may  be 
compared  by  examining  Figures  7,  8,  18,  and  19.  The  data  presented  in 
these  figures  are  for  the  pitch  axis  located  at  the  quarter  chord  and  for 
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oscillations  about  small  mean  angles  of  attack.  Theoretical  results  are 
presented  in  all  figure*  so  that  the  data  may  be  compared  by  noting  the 
agreement  with  theory.  The  results  of  the  two  models  show  excellent 
agreement  except  in  the  case  of  pitching  moment  phase  angle,  as  noted 
previously. 

Instantaneous  normal  force  and  pitching  moments  are  presented  for  the  two 
models  in  Figures  9,  10,  and  29.  The  instantaneous  curves  for  the  two 
models  have  similar  shape  and  indicate  the  same  trends  with  increasing 
oscillating  frequency.  The  differences  in  shape  of  the  individual  curves 
obtained  for  a  given  oscillating  frequency  are  to  a  large  measure  attributed 
to  the  differences  in  the  stall  characteristics  of  the  two  models  and  the 
relationship  of  the  mean  angle  of  attack  to  the  steady  state  stall  angle  of 
attack.  The  minimum  pitching  moment  values  for  the  pressure  model  and 
the  force  model  of  -0.  33  and  -0.  36,  respectively,  compare  exceedingly 
well.  Also,  from  Figures  9  and  29,  for  which  the  mean  angle  is  close  to 
the  stall  angle,  the  maximum  Cjj  values  obtained  for  both  models  of  approxi¬ 
mately  1. 8  agree  quite  well. 


Comparison  With  Previous  Results 

Some  previous  investigations  are  reported  in  Appendix  I  for  airfoils  oscilla¬ 
ting  at  low  mean  angles  of  attack.  Once  again  the  data  are  compared  with 
theory  and  can  be  compared  with  the  results  of  this  investigation  using 
theoretical  values  as  a  guide. 

Instantaneous  normal  force  and  pitching  moment  coefficients  for  airfoils 
oscillating  in  the  stall  region  are  presented  by  HalfmanJ,  Carta6,  and 
Liiva,  Davenport,  Gray  and  Walton7.  Halfman  presents  data  for  three 
12-percent-thickness  airfoils  which  he  refers  to  as  sharp,  blunt,  and  inter¬ 
mediate.  The  intermediate  airfoil  is  very  similar  to  the  NACA  0012  airfoil. 
Carta  tested  a  0012  airfoil.  Liiva  et  al  tested  a  Vertol  23010-1.  58  airfoil 
and  a  0012  airfoil. 

Data  presented  by  Halfman  are  for  a  pitch  axis  at  37  percent  chord.  Only 
instantaneous  pitching  moments  are  presented,  and  only  a  portion  of  this 
is  for  the  intermediate  airfoil.  The  data  presented  for  comparable  condi¬ 
tions  show  essentially  the  same  type  loops  and  indicate  the  large  drop  in 
pitching  moment  as  reported  herein.  The  two  investigations  were  conducted 
at  the  same  value  of  Reynolds  number.  Liiva  et  al  present  both  instanta¬ 
neous  normal  force  and  pitching  moment  coefficients  for  an  NACA  0012  air¬ 
foil  oscillating  about  the  quarter  chord  at  a  Mach  number  of  0.4.  With  the 
exception  of  Mach  number,  oscillating  conditions  are  nearly  the  same  as 
those  reported  herein  for  the  25  percent  pitch-axis  models.  The  higher 
Mach  number  changes  the  steady  state  data  somewhat,  but  the  shape  of  the 
Cn  versus  angle-of-attack  curve  at  the  stall  point  is  very  similar  to  that  of 


24 


the  force  model  of  this  report.  The  curves  show  good  agreement  in  shape, 
magnitude, and  trends.  Data  of  this  repor+  are  for  a  mean  angle  closer  to 
the  stall  angle  than  that  reported  by  Liiva.  Therefore,  the  model  approaches 
the  unstalled  condition  at  a  lower  value  of  frequency  parameter,  causing  the 
curves  of  this  report  to  agree  with  those  at  a  somewhat  higher  value  of  k  in 
Liiva's  report.  The  values  of  Cj^  maximum  agree  very  well  (1.7  as  com¬ 
pared  to  1.  8  of  this  report).  The  pitching  moment  drop  shows  excellent 
agreement  in  magnitude  (-0.  33  as  compared  to  -0.  36  of  this  report). 

Halfman®  presents  some  drag  data  for  an  NACA  0012  airfoil  oscillating  in 
pitch  about  the  37  percent  chord  at  a  mean  angle  of  zero  degrees.  He  pre¬ 
sents  an  average  drag -amplitude  coefficient  as  a  function  of  k.  This  drag- 
amplitude  coefficient  increases  with  increasing  k.  This  trend  is  contradic¬ 
tory  to  the  results  of  this  investigation,  which  indicate  that  the  mean  value 
increases  with  increasing  k  but  the  oscillatory  magnitude  remains  essentially 
constant  (a  -  0  degrees)  or  decreases  with  increasing  k  (a  =  6  degrees). 

There  are  no  other  oscillatory  drag  data  available  to  support  either  of  these 
investigations. 
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CONCLUSIONS 


The  literature  survey  revealed  considerable  low-speed  experimental  data 
on  oscillating  two-dimensional  airfoils.  Very  little  of  this  data  can  be 
directly  correlated  due  to  the  number  of  parameters  involved  (oscillating 
frequency,  oscillating  amplitude,  mean  angle  of  attack,  test  Reynolds 
number,  airfoil  profile  and  pitch-axis  location)  and  the  choice  of  data 
recorded  by  the  investigators.  In  cases  where  direct  correlation  is  possi¬ 
ble,  scatter  of  the  data  restricts  the  possibility  of  any  definite  conclusions. 

The  agreement  between  experimental  data  of  this  investigation,  obtained 
from  the  pressure  measurements  and  the  direct  force  measurements,  and 
the  theoretical  data  at  low  mean  angles  of  attack  indicates  that  both  methods 
can  produce  satisfactory  results  over  the  frequency  range  tested.  How¬ 
ever,  it  is  felt  that  the  present  force  balance  was  operating  very  near  its 
useful  limit  at  the  higher  oscillating  frequencies.  Both  methods  have 
their  advantages  and  disadvantages.  The  advantages  of  each  are  as  follows: 

Pressure  measurements: 

1.  Results  are  not  affected  by  inertia  loads. 

2.  Model  supports  can  be  rigid  enough  that  natural  frequencies 
present  no  problems. 

3.  Instantaneous  pressure  distributions  may  be  obtained. 

Force  measurements: 

1.  Requires  less  instrumentation. 

2.  Less  data  processing. 

3.  Possible  to  obtain  drag  information.  (Oscillating  drag  trends 
were  obtained  in  this  investigation, but  actual  drag  magnitudes 
are  questionable  due  to  lack  of  sensitivity  and  low  natural 
frequency  of  the  drag  balance.) 

The  disadvantages  are: 

Pressure  measurements: 

1.  Drag  data  are  not  available. 

2.  Requires  considerable  instrumentation  for  a  reasonable 
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number  of  pressure  orifices. 

3.  Requires  considerable  data  processing. 

4.  Difficult  to  install  sufficient  number  of  pressure  transducers 
in  model. 

Force  model: 

1.  Inertia  loads  must  be  transmitted  through  the  force  balance. 
This  imposes  problems  for  balance  design. 

2.  Flexibility  required  for  balance  sensitivity  results  in  low 
natural  frequencies. 

3.  Inertia  loads  must  be  eliminated  from  measurements  in  some 
way. 

The  fact  that  the  pressure  data,  which  were  measured  along  the  centerline 
of  the  model  where  the  gaps  between  the  model  and  the  tunnel  walls  had 
minimal  effect  upon  the  results,  agree  well  with  the  force  data  indicates 
that  the  small  gaps  employed  during  the  test  had  little  influence  upon  the 
data. 

From  this  investigation,  it  may  be  concluded  that: 

1.  For  values  of  instantaneous  angle  of  attack  not  exceeding  the 
steady  state  stall  value,  theoretical  predictions  show  good  agree¬ 
ment  with  test  results. 

2.  When  oscillating  at  high  frequency  about  mean  angles  above  the 
steady  state  stall  angle,  values  of  Cfl  much  greater  than  the 
steady  state  maximum  value  may  be  obtained.  A  value  of  2.  29 
was  obtained  during  this  investigation.  This  is  approximately  80 
percent  greater  than  the  steady  state  value.  It  may  be  possible  to 
achieve  considerably  higher  values  with  appropriate  values  of  a, 
and  k. 

3.  For  the  models  with  the  pitch  axis  at  the  quarter  chord,  sharp 
decreases  in  pitching  moment  were  experienced  at  some  fre¬ 
quencies  when  the  instantaneous  angle  of  attack  exceeded  the 
steady  state  stall  values.  The  minimum  value  noted  during  this 
investigation  was  -0,36.  The  decreases  in  pitching  moment  were 
less  severe  as  the  pitch  axis  was  moved  toward  the  midchord. 

4.  Mean  drag  values  increase  rapidly  with  increasing  oscillating  fre¬ 
quency. 
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5.  Drag  results  of  this  investigation  indicate  a  need  for  better 
drag  studies. 
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RECOMMENDATIONS 


It  is  recommended  that  the  possibility  of  increasing  the  sensitivity  of  the 
drag  balances,  by  use  of  semiconductor  strain  gages,  and  increasing  the 
natural  frequency  of  the  balance  be  investigated  to  obtain  more  accurate 
oscillating  drag  data. 

The  unusual  pressure  distributions  noted  on  the  lower  surface  leading  edge 
of  the  model  under  some  oscillating  conditions  indicate  a  possible  vortex 
formation.  It  would  be  desirable  to  do  some  flow  visualization  studies  to 
determine  the  nature  of  the  flow  in  this  region. 
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FIGT  RE  1.  Oscillating  Drive  Mechanism  Mounted  Under  the  Test  Sec¬ 
tion  of  the  1.  5-Foot  by  3.  87  5-Foot  Wind  Tunnel. 
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FIGURE  2.  View  With  Tunnel  Wall  Removed  of  Force  Model  (Axis  at  50%  Chord)  Installed 
in  Test  Section. 


FIGURE  3.  Details  of  Pressure  Transducer  Installation. 
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FIGURE  5.  Strain  Gage  Electrical  Connections. 
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FIGURE  6.  Static  Cj^  and  Characteristics. 
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FIGURE  7 
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Variation  of  Normal  Fores  Coefficient  With  Reduced  Fre¬ 
quency  for  Pressure  Model  Oscillating  in  Pitch,  Pitch  Axis 
•t  25%  Chord. 
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FIGURE  8.  Variation  of  Pitching  Moment  Coefficient  With  Reduced  Fre¬ 
quency  for  Pressure  Model  Oscillating  in  Pitch,  Pitch  Axis 
at  25%  Chord. 
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FIGURE  9.  Continued. 
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FIGURE  0.  Continued. 
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FIGURE  10.  Continued. 
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Instantaneous  Pressure  Distributions,  a  =  10.  05°,  5=5. 
=  6.0°,  k  =  .  032. 


4fi 


PRESSURE  COEFFICIENT  (C 


&a  =  6.  17°,  k  =  .  268. 


47 


PRESSURE  COE 


FIGURE  13.  Continued. 
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FIGURE  16.  Instantaneous  Pressure  Distributions,  o  *  18.00*, 
Ao  s  t>17*.  k  «  .*66. 
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FIGURE  18.  Variation  of  Normal  Force  Coefficient  With  Reduced  Fre¬ 
quency  for  Force  Model  Oscillating  in  Pitch.  Pitch  Axis 
at  25%  Chord,  &a  *  6.08*. 
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FIGURE  19.  Variation  of  Pitching  Moment  Coefficient  With  Reduced 
Frequency  for  Force  Model  Oscillating  in  Pitch,  Pitch 
Axis  at  2  5%  Chord,  &a  =  6.  08°. 
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FIGURE  20.  Variation  of  Normal  Force  Coefficient  With  Reduced  Fre¬ 
quency  for  Force  Model  Oscillating  in  Pitch,  Pitch  Axis 
at  37%  Chord,  &a  =  6.08°. 
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FIGURE  22.  Variation  of  Normal  Fore*  Coefficient  Uith  Rcojced 
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Variation  of  l  itching  Moment  Coefficient  With  Reduced 
Frequency  for  Force  Model  Oscillating  in  Pitch,  Pitch 
Axis  at  50%  Chord,  A  a  =  6.08°. 
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FIGURE  24.  Variation  of  Drag  Coefficient  With  Reduced  Frequency, 
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FIGURE  25.  Variation  of  Drag  Coefficient  With  Reduced  Frequency, 
a  ss  6°. 


66 


0.05 


0.04 


0.03 


0.02 


0.01 


0 


ANGLE  OF  ATTACK  (cr)-DEGREES 

FIGURE  2  7.  Dynamic  CN  and  CM  for  50%  Pitch  Axis  Model  Oscillating 
at  High  Frequency  About  a  =  6.  22°,  k  =  .  297,  &a  =  6.  38°. 
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FIGURE  29.  Continued 
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3 


PITCHING  MOMENT  COEFFICIENT  (C,,)  NORMAL  FORCE  COEFFICIENT  (C 


PITCHING  MOMENT  COEFFICIENT  (CM>  NORMAL  FORCE  COEFFICIENT  (C 
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FIGURE  30.  Continued 
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APPENDIX  I 

CORRELATION  OF  OSCILLATING  AIRFOIL  DATA 


In  attempting  to  compare  experimental  data  from  various  sources,  two 
problem  areas  exist.  The  first  of  these  is  the  large  number  of  parameters 
involved:  oscillating  frequency,  oscillating  amplitude,  test  Reynolds  num¬ 
ber,  airfoil  profile,  and  pitch-axis  location.  The  second  problem  area  is 
in  the  type  of  data  reported.  That  is,  some  investigators  present  only  lift 
information,  whereas  others  present  only  moment  data.  Oddly  enough, 
early  work  in  this  country  was  concerned  only  with  lift  in  pure  pitch;  where¬ 
as  the  British  investigators  measured  only  moments  in  pure  pitch.  A 
study  of  Table  II,  a  tabular  summary  of  experimental  work,  will  point  out 
the  extent  of  the  problems  encountered  in  correlating  or  comparing  experi¬ 
mental  results. 

Theoretical  values  of  the  aerodynamic  coefficients  provide  one  basis  for 
comparison  of  the  experimental  results.  Linearised  theory  applies  to 
very  thin  wings  undergoing  small  oscillations  in  an  inviscid  fluid.  With 
these  limitations  in  mind,  comparison  may  be  made  between  theory  and  low- 
speed  investigations  of  airfoils  oscillating  about  low  mean  angles  of  attack. 
For  purposes  of  this  report,  the  theory  of  Theodorsen^  has  been  used  as  a 
guide. 

Experimental  work  on  oscillating  airfoils  is  complicated  by  the  difficulty 
in  obtaining  pure  harmonic  motions  and  the  problems  associated  with 
measuring  the  forces  and  moments.  In  addition  to  these  obstacles,  calcu¬ 
lation  of  tunnel  wall  corrections  is  a  formidable  task,  so  most  data  are 
presented  uncorrected.  The  experimental  difficulties  coupled  with  the 
theoretical  approximations  result  in  a  standard,  of  what  is  considered  to 
be  good  agreement,  being  set  rather  low.  Results  presented  in  the  investi¬ 
gations  of  this  survey  for  airfoils  operating  at  small  mean  angles  of  inci¬ 
dence  (such  that  the  stall  angle  is  not  exceeded)  show  qualitative  agreement 
with  the  linearized  theory. 

O 

Halfman  presents  an  excellent  survey  of  the  first  six  references  of  Table 
II.  In  order  to  compare  the  data,  he  had  to  modify  them  to  agree  in  form 
and  correct  for  the  differences  in  pitch-axis  location  used  by  the  various 
investigators.  In  correcting  for  the  pitch-axis  location.  Halfman  had  to 
resort  to  theoretical  considerations  since  there  were  insufficient  experi¬ 
mental  data  to  perform  the  required  transfers.  After  applying  the 
corrections,  there  was  still  a  large  difference  in  Reynolds  number.  Lift 
and  moment  data  were  then  plotted  against  Reynolds  number  for  various 
reduced  frequencies  to  provide  a  comparison.  These  plots  from  Halfman's 
report  have  been  reproduced  here  as  Figures  32  and  33.  In  Figure  32  the 
curves  labeled  (S)  are  from  the  work  of  Reid  and  Vencenti*!  and  Reid^, 
those  labeled  (M)  from  the  work  of  Halfman^,  and  those  labeled  (T)  from 
the  theory  of  Theodorsen^.  Considering  the  corrections  required  to 
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modil>  the  data  and  the  fact  that  the  airfoils  are  slightly  different  profiles, 
figure  32  shows  remarkable  agreement  in  the  trends.  In  Figure  33,  the 
curves  Labeled  (B|)  are  from  the  data  of  Bratt  and  Scruton®,  those  labeled 
(B2)  from  Bratt  and  Wight those  labeled  (M)  from  Halfman**,  and  those 
labeled  (T)  from  theory.  In  the  moment  phase  anele  portion  of  Figure  33, 
it  may  he  noted  that  the  data  of  Reid  and  Vencenti11  and  Reid  ^  show  re¬ 
markable  agreement,  where,  s  tnose  of  Bratt  and  Scruton®  do  not  agree  as 
well.  The  moment  magnitude  data  are  not  as  consistent  as  the  moment 
phase  angle  data.  An  interesting  point  in  the  moment  magnitude  data  is 
that  a  sharp  increase  in  magnitude  with  Reynolds  number  occurs  at  the 
higher  reduced  frequencies  for  all  three  sets  of  data  presented.  No  expla¬ 
nation  of  this  phenomenon  is  offered.  Again,  it  should  be  borne  in  mind 
that  the  data  from  the  three  sources  compared  in  Figure  33  are  not  for  the 
same  airfoil  profiles. 

In  an  attempt  to  determine  the  effects  of  airfoil  profile,  lift  and  moment 
data  from  Halfman**  and  Halfman  et  ar  have  been  compared  in  Figures  34 
and  35  respectively.  Data  presented  are  for  four  different  profiles  of  the 
same  thickness  ratio  obtained  from  the  same  test  facility.  Whereas 
qualitative  agreement  with  theory  may  be  noted,  the  scatter  of  the  data 
precludes  any  possibility  of  determining  the  effects  of  airfoil  profile. 

Along  this  same  line,  data  from  Wyss  and  Monfort17  are  presented  in 
Figures  36  and  37.  The  data  presented  are  for  five  airfoils  with  thickness 
and  chordwise  location  of  maximum  thickness  varied  systematically. 
Although  the  data  presented  are  for  a  Mach  number  0.491,  the  compressi¬ 
bility  effects  should  be  quite  small  since  the  mean  angle  of  incidence  (2°) 
and  the  oscillating  amplitude  (1*)  are  pmall.  Once  again,  there  is  qualita¬ 
tive  agreement  with  theory  but  too  much  scatter  in  the  data  to  arrive  at 
any  conclusions  regarding  the  effects  of  orofile. 

No  other  attempt  has  been  made  to  compare  the  data  from  the  surveys 
listed  in  Table  11  because  of  the  lack  in  conformity  of  the  parameters  in¬ 
volved  and  the  fact  that  each  investigation  is  compared  with  theory  in  the 
individual  reports. 

Three  references (5,  6,  and  9)lisled  in  Table  II  present  data  for  airfoils 
oscillating  at  high  angles  of  attack  such  that  the  airfoil  is  operating  in  the 
stalled  region  part  of  the  time.  It  may  be  noted  from  the  table  that  the 
pitch -axis  location  is  different  for  each  investigation.  Since  there  is  no 
theoretical  correction  which  can  be  applied  to  account  for  the  different 
pitch-axis  location  for  airfoils  oscillating  at  large  angles  of  attack,  the 
data  cannot  be  compared  directly.  While  References  5  and  9  present 
moment  data,  besides  the  difference  in  pitch-axis  location,  they  are  for 
different  airfoil  shapes  and  different  Reynolds  numbers.  Data  presented 
from  the  investigations  reported  in  References  5  and  6  were  for  the  same 
airfoil  shape.  Whereas  normal  force  is  presented  for  instantaneous  angles 
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of  attack  in  Reference  6,  corresponding  data  arc  not  presented  in  Reference 
5.  So  here  again,  there  is  no  chance  for  direct  comparison. 

In  addition  to  the  tabular  summary  (Table  II)  of  two-dimensional,  low- 
speed,  experimental  oscillating  airfoil  investigations,  tabular  summaries 
are  presented  in  Tables  ’Il-V  for  experimental  investigations  of  two- 
dimensional  oscillating  airfoils  in  compressible  flow  and  finite  wings  oscil¬ 
lating  in  low  speed  and  compressible  flow.  Tables  111  -  V  are  not  intended 
to  represent  a  complete  survey.  They  are  a  summary  of  investigations 
that  were  brought  to  attention  during  the  search  for  low-speed  two-dimen¬ 
sional  investigations  and  are  included  in  this  report  for  convenience.  No 
attempt  is  made  to  correlate  these  summaries. 

In  conclusion,  while  there  is  considerable  low  -speed  experimental  data  on 
oscillating  two-dimensional  airfoils,  there  is  very  little  that  can  be  directly 
correlated.  This  is  due  to  the  number  of  parameters  involved  and  the 
choice  of  data  recorded  by  the  investigators.  In  cases  where  direct  corre¬ 
lation  is  possible,  scatter  of  the  data  restricts  the  possibility  of  any 
definite  conclusions.  Whereas  qualitative  agreement  with  linearized  theory 
is  indicated  in  most  cases,  quantitative  agreement  is  hampered  by  the 
theoretical  assumptions  and  the  difficulty  of  obtaining  accurate  experimen¬ 
tal  data. 
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TABLE  II.  OSCILLATING  AIRFOIL,  SUMMARY  OF 


Her. 

No. 

Air  foil 

Sect  ion 

Chord 

In. 

Pitch 

Axis  %  C 

0 

Deg. 

Act 

Deg. 

Trans. 

Amp.  In. 

9 

15%  .Toukowski 

9 

.  09 -. 28 

50 

0-20 

2.  1-10 

N/A 

10 

18%  Symmetrical 

5.  188 

. 06-. 3 

25 

0 

N/A 

1 1 

NACA  0015 

15 

.  I-?-.  30 

40 

0,  -5 

2. 5-7. 5 

N/A  4 

12 

NACA  0015 

10  &  15 

.  137-1. 028 

30  &  40 

0-10 

1-5 

N/A 

13 

15%  Joukowski, 

15%  EC  1550,  15% 
Elliptical,  15% 
Hollow -Ground 

9 

.  i42-  . 283 

33. 3  &  50 

0-18 

2-6 

N/A 

8 

NACA  0012 

12 

.7  -  1.0 

37 

0,  6.  1 

5.  19,6.7  &  13.  5 

1.0,  1.37,2.0 

5 

12%  Symmetrical 
Sharp  Intermediate 
Blunt 

12 

.7  -1.0 

37 

0-22 

6.08 

0.9 

14 

7.  3%  Symmetrical 

11.7 

.08  -1.0 

Varied 

0 

.95-1.9 

.63 

6 

NACA  0012 

24.  5 

3.  0  -  6. 0 

25 

0-33 

4-8 

N/A 

15 

NACA  0018 

4 

0.2-1. 2 

25 

0 

.5 

N/A 

36 

NACA  632A613 

8 

18.75  & 

32.  1 

9.  0 

. 78-. 9fc 

N/A 
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[ARY  OF  LOW -SPEED  EXPERIMENTAL  INVESTIGATIONS 


Data  Presented 


ins. 

3.  In. 

Freq. 

CPS 

k 

Instrumentation 

Lift 

N. 

Force 

Drag 

Moment 

Phase 

Angle 

Remarks 

(A 

0-11 

0-.  75 

Magneto-Striction 
Stress  Indicator 

X 

Virtual  Mass 
Effect  of  Air  is 
Cancelled  Out  of 
Data 

f A 

17 

.  07-.7 

Spring  Deflection 

X 

t A 

4.24  &  8.82 

.2  -.95 

Dynamometer 

X 

X 

>A 

6.6-15 

.2-2.0 

Dynamometer 

X 

X 

>A 

0-11 

0-  .  785 

Magneto-Striction 
Stress  Indicator 

X 

.37,2.0 

0-17 

.05-.  46 

Strain  Gage 

X 

X 

X 

X 

Combined  Pitch 
and  Translation 

9 

2-18 

.05-.  55 

Strain  Gage 

X 

X 

X 

X 

63 

4.  17  &  8.33 

.  08-1.0 

Electro-Dynamic 

Pickup 

X 

X 

Combined  Pitch 
and  Translation 
in  Kept.  F.  102 

'A 

0-16 

.06-  .3 

Pressure  Trans¬ 
ducers 

X 

X 

'A 

15 

.04-2.4 

Piezo-Electric 

Gages 

X 

X 

X 

Water  Tunnel 

rA 

18.61-19.5 

.746- 

2.66 

Strain  Gage 

X 

X 

i 


TABLK  III.  OSCILLATING  AIRFOIL,  SUMMARY  OF  COMPRESSIBLE  FLOW  1 


Ref. 

A  irfoil 

( 'lion! 

Mach 

KN/I0(l 

Bitch 

a 

Ac 

Trans. 

Freq. 

N». 

Sect  ion 

In. 

No. 

Range 

Axis  %  C 

Deg. 

Deg. 

Amp.  In. 

CPS 

k 

Instrumentation 

16 

Bi-Convex 

2 

0.  4-0. 9 
l.  275 

.  4-.  8 

50 

0 

2-4 

N/A 

8-27 

0-.  015 

Photocell 

17 

NACA 

24 

.  5-  .  9 

5. 0-8. 0 

25 

o,  2 

1 

N/A 

4-40 

.  025- 

Pressure 

65A012 

65A008 

2-008 

877A008 

0.  45 

Transducers 

65A004 

4 

NACA 

24 

.2-. 86 

3. 0-8. 0 

25 

4-  10 

1 

N/A 

4-40 

.03- 

Pressure 

65A012 

65A008 

2-008 

877A008 

1.  12 

Transducers 

65A004 

18 

NACA 

24 

.5-.  9 

5. 0-8. 0 

25 

2 

1 

N/A 

4-40 

.025- 

Pressure 

65A012 

65A008 

2-008 

877A008 

0.  45 

Transducers 

19 

NACA 

12 

. 35-. 70 

5.  0 

50 

0 

1.2, 

N/A 

0-60 

0-.  8 

Pressure 

65-010 

2.4 

Transducers 

20 

NACA 

12 

.  35-.  70 

5.  3 

50 

0-16 

1.2 

N/A 

10-35 

0-.64 

Pressui  e 

65A010 

Transd  era 

21 

NACA 

65-010 

8 

. 35-. 78 

1.  0- 

25 

0 

2.  3 

N/A 

.  17- 

Strain  Gages 

5.  5 

.35 

Wattmeter 

22 

Blunt-Nosed 

2.  5 

1.  75- 

0.  9 

Varied 

0 

N/A 

20 

.  12 

Light  Beams 

Wedge 

2.  47 

28 

Double 

3 

8.  8 

1.  0 

50 

0 

3 

N/A 

125 

.018 

Strain  Gage 

Wedge 

24 

Double 

1.  37- 

Varied 

0 

1 

N/A 

Wedge 

2.  43 

1— 

Single  Wedge 

AIRFOIL.  SUMMARY  OF  COMPRESSIBLE  FLOW  EXPERIMENTAL  INVESTIGATIONS 


Data  Presented 

Trans. 
Amp.  In. 

Freq. 

CPS 

k 

Instrumentation 

Lift 

N. 

Force  Drag  Moment 

Phase 

Angle 

Damping 

Remarks 

N/A 

8-27 

0-. 015 

Photocell 

X 

Measured  Model  Dis¬ 
placement 

N/A 

4-40 

.  025- 
0.45 

Pressure 

Transducers 

X 

X 

t 

X 

N/A 

4-40 

.  03- 
1.  12 

Pressure 

Transducers 

X 

X 

X 

N/A 

4-40 

.  025- 
0.45 

Pressure 

Transducers 

X 

X 

X 

Effect  of  Spoilers  on 
Oscillating  Airfoils 

N/A 

0-60 

0-.  8 

Pressure 

Transducers 

X 

X 

X 

N/A 

10-35 

CO 

1 

o 

Pressure 

Transducers 

X  X 

X 

N/A 

.  17- 
.  35 

Strain  Gages 
Wattmeter 

X 

X 

X 

N/A 

20 

.  12 

Light  Beams 

X 

Measured  Model  Dis¬ 
placement 

N/A 

125 

.  018 

Strain  Gage 

N/A 

X 

Ref.  Airfoil  Chord 

No.  Section  In.  AR  TK  c/4 


15%  Thick 

10 

1.0 

0 

6 

1.0 

0 

4 

1.0 

0 

2 

1.  0 

0 

1 

1.0 

0 

15%  Joukow- 
ski 

9 

2.7 

1.0 

0 

15%  Joukow- 
ski,  EC  1550 

9 

1.0 

1.0 

0 

2.  0 

1.0 

0 

3.  0 

1.0 

0 

3.  5 

1.0 

0 

4.4 

1.  0 

0 

NACA  0020 

3.  75 

3 

1.0 

0 

3.  75 

4 

1.0 

0 

3.  75 

5 

1.0 

0 

3.  75 

3 

1.0 

45 

3.  75 

5 

1.0 

45 

8%  Gothic 

.  75 

5%  Gothic 

.  75 

RAE  102 

1.2 

14 

RAE  102 

24 

1.6 

0 

RAE  102 

21 

1.  32 

46 

RAE  102,  10% 

13.  4 

2.  97 

.  14 

36.  & 

RAE  101,  6% 

21.3 

3.0 

1.0 

60 

EQ  10,  40.  10% 

20.  7 

4.4 

.31 

40 

RAE  101 

20 

3.  3 

1.  0 

0 

MEASURED  FROM  WING  APEX 


T  ABLE  IV.  OSCILLATING  WING,  SUMMARY  OF  LOW-SPE 
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Pitch  Axis 

a 

Ac 
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%  C 

Deg. 
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0.  9 

50 

24.  3 

.  19- 

50 

0.6 
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.  08 
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50 
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6 
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18 

10 

27 
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* 
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60 
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28.  8 

.75- 

# 

1 

.5 
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0-15 
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N/A 
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33.6  1. 

2 
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2-4.  7 

38.5  1. 

7 

1. 112* 

0-15 

1.  15- 

.883* 

5.  15 

40.2 

.7-2.2 

.328* 

0 

1.65 

1.0- 

.055 

2.0 

64 

.288 

0 

1.70 

97.  5 

.258 

0 

1.70 

.  572* 

.  4-1.  5  0 

100 


33.  5 


0 


2. 0  N/A 


ING,  SUMMARY  OF  LOW -SPEED  EXPERIMENTAL  INVESTIGATIONS 


Data  Presented 

Axis  a 

Ac 

Trans. 

Freq. 

Instru- 

Phase 

Deg. 

Deg. 

Amp.  In. 

CPS 

k 

mentation 

Lift 

Moment 

Angle 

Damping 

Remarks 

0.  35 

X 

X 

X 

Pitch  and 

Heave 

0.6 

5-10 

N/A 

0-10. 3 

0- .  506  Magneto 

X 

Striction 
Stress  In¬ 
dicator 

0.8. 

6 

N/A 

0-11 

0-.  78 

Magneto 

X 

10 

Striction 
Stress  In¬ 
dicator 

N/A 

N/A 

0.8 

.025- 

Light  Beam 

X 

Measured  Model 

.25 

Displacement 

0-20 

N/A 

r- 

i 

o 

X 

X 

0-15 

.9-4.4 

N/A 

.  03- 

.375 

Light  Beams  X 

X 

Clipped  Delta 

0-15 

2-4.  7 

.  03- 

Wing  6%,Delta 

.375 

Wing  6% 

0-15 

1.  15- 

.  03- 

Arrowhead  wing 

5.  15 

.  375 

10% 

0 

1.65 

1.0- 

0-2.97 

o 

1 

o 

00 

Strain  Gages  X 

X 

X 

Only  90*  Delta 

2.0 

Wing  Heaved 

0 

1.70 

0-2.97 

0-.  13 

0 

1.70 

0-2.97 

0-.  13 

Measurements 

Made  Using  Free 
and  Forced  Os- 

0 

2.0 

N/A 

.2-. 65 

Force  Trans- 
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ducer 

TABLE  V.  OSCILLATING  WING,  SUMMARY  OF  CO 


Wing 

Charac 

■t  eristics 

Ref. 

Airfoil 

( ’hord 

Semi  - 

Mach 
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a 

Aa 
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No. 
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In. 
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%  C 

Deg. 

Deg. 

Amp.  I 

31 
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5 

5.  4 

1.0 

0 
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.  049- 

25 

0-24 

3 
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5.  04 

1. 0 

15 
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30 
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50 
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45 
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75 
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65A004 

8 
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0 

24 

8 
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0 

21 

8 
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18 

8 
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15 
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0 
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50 
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50 

0 
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40 
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0 
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0 

.8- 

N/A 
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* 
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* 
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3 
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IMARY  OF  COMPRESSIBLE  FLOW  EXPERIMENTAL  INVESTIGATIONS 


Data  Presented 


A  a 
Deg. 

Trans. 
Amp.  In. 

Freq. 

CPS 

Instrumen- 
k  tation 

Phase 

Lift  Moment  Angle 

Damping 
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3 

N/A 

35-175 

.15-1.3  Strain  Gages 

X 

Flutter  Tests 

N/A 
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X 

X 
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Wing  With  Tip  Tank 
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2.  11 

N/A 
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.  15-1.  32 
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X 

X 

X 

1.5 

N/A 
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.  008- 
.269 

Pressure 

Trans. 

X 

X 

X 

.8- 
.  6 

N/A 

35-70 
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.  187 

Dampo- 

meter 

X 
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REDUCED  FREQUENCY  00 

FIGURE  34.  Airfoil  Profile  Effect  -  Lift  in  Pure  Pitch,  a  -  0#, 
&a  ■  6.08*.  Pitch  Axis  «  37%  Chord. 


113 


PHASE  ANCLE  (*  )  -  DEGREES  MAGNITUDE 
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FIGURE  37.  Effect  of  Thickness  Distribution  -  Lift  in  Pure  Pitch, 
M  =  .491,  a  =  2°.  Pitch  Axis  =  25%  Chord. 
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APPENDIX  II 

EFFECT  OF  TUBING  ON  THE  REMOTE 
READING  OF  OSCILLATING  PRESSURES* 


INTRODUCTION 

In  order  to  support  the  investigation  of  the  pressure  distribution  over  an 
oscillating  airfoil,  it  is  desirable  to  have  a  knowledge  of  the  effect  of  tubing 
size  and  diameter  upon  the  response  of  a  remotely  located  transducer  at 
various  frequencies  of  airfoil  oscillation.  There  are  two  aspects  to  the 
response  --  attenuation  of  the  pressure  and  pressure  lag.  For  this  investi¬ 
gation,  the  term  pressure  attenuation  is  defined  as  the  ratio  of  the  remote 
pressure  to  the  local  or  reference  pressure. 

It  was  hoped  that,  in  order  to  keep  the  amount  of  recording  equipment  to  a 
minimum,  lengths  of  tubing  could  be  led  from  the  pressure  orifices  on  the 
wind-tunnel  model  to  a  scanivalve  and  single  transducer  located  outside  the 
wind  tunnel. 

During  the  course  of  the  investigation,  it  was  decided  that  it  would  be  use¬ 
ful  to  have  an  idea  of  the  error  encountered  when  transducers  are  mounted 
inside  the  model  and  connected  to  the  orifices  by  short  lengths  of  tubing. 

EXPERIMENTAL  APPARATUS  AND  PROCEDURE 

The  apparatus  consisted  of  a  sinusoidally  oscillating  piston  which  was 
driven  by  a  Scotch-yoke  mechanism.  Two  differential  pressure  transducers 
(Statham  Model  No.  PM  131TC  ±2.  5-350)  were  used,  one  whose  diaphragm 
was  mounted  flush  with  the  end  of  the  piston  cylinder  and  another  that  was 
connected  to  the  cylinder  through  lengths  of  flexible  tubing.  The  signals  of 
both  pressure  transducers  were  first  fed  into  a  multichannel  carrier 
amplifier  and  thence  into  a  recording  oscillograph.  Calibration  of  the  re¬ 
cording  system  was  accomplished  by  comparing  readings  from  a  water - 
filled  manometer  with  stepped  oscillograph  traces. 

The  experimental  investigation  consisted  of  three  phases.  The  first  phase 
examined  the  effects  of  tubing  length  and  oscillation  frequency  --  the  re¬ 
mote  transducer  being  mounted  in  a  scanivalve.  The  second  phase  examined 
the  effects  of  tubing  length  and  frequency  of  oscillation  amplitude  -  the  re¬ 
mote  transducer  being  mounted  in  a  special  holder.  The  third  phase 
examined  the  effects  of  tube  diameter,  frequency  of  oscillation,  oscillation 


*Dagold,  Reuben  G. ,  RESEARCH  MEMORANDUM  1967-4,  Wind  Tunnel 
Operations  Department,  University  of  Maryland,  College  Park,  Maryland, 
1967. 
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amplitude,  and  tubing  connection  position  (i.  e. ,  to  the  flush  diaphragm  face 
of  the  transducer  and  to  the  plug  side  of  the  diaphragm)  with  a  short  length 
of  tubing  --  the  remote  transducer  again  being  mounted  in  the  special  holder. 

The  repeatability  of  the  data  between  runs  using  the  scanivalve  (first  phase) 
was  poor  (see  Figures  38  and  39).  A  very  flexible  rubber  gasket  is  used 
to  seal  the  inlet  tubes  of  the  scanivalve.  The  gasket  can  partially  restrict 
the  opening;  the  extent  of  the  restriction  changes  each  time  that  the  seal  is 
tightened  and  also  varies  from  port  to  port.  Thus,  it  was  surmised  that 
the  poor  repeatability  might  be  attributed  to  the  scanivalve  design. 

A  special  holder  was  made  that  avoided  the  varying  restriction  problem  of 
the  scanivalve.  Use  of  the  holder  during  the  second  phase  produced  good 
repeatability  of  data  (see  Figures  40  and  41). 

The  entire  investigation  was  conducted  at  a  mean  pressure  level  of  one 
atmosphere,  that  being  the  pressure  in  the  cylinder  where  the  piston  is  at 
midstroKe.  It  was  at  that  piston  position  and  pressure  that  phase-lag  data 
were  obtained.  The  phase  lag  is  defined  here  as  the  amount,  in  degrees, 
that  the  remotely  read  pressure  lags  the  pressure  recorded  in  the  cylinder. 
By  amplitude,  one-half  the  difference  between  the  maximum  and  minimum 
pressures  is  meant. 

DISCUSSION  OF  RESULTS 


It  was  noted  that  phase  lag  varies  during  a  cycle  and  that  although  the  pis¬ 
ton  motion  was  sinusoidal,  the  pressure  did  not  vary  sinusoidally.  Figure 
42  provides  an  example  of  lag  variation  over  the  cycle  (the  curves  have 
been  normalized).  It  is  evident  here  that,  for  the  remote  pressure  reading, 
the  time  for  the  pressure  to  decrease  from  the  maximum  to  the  minimum  is 
greater  than  the  time  for  the  pressure  to  increase  from  the  minimum  to  the 
maximum.  This  is  substantiated  by  Larcombe  and  Peto^,  who  state  that 
for  a  constant  absolute  value  of  difference  between  the  initial  pressure  in 
the  tube  and  the  applied  pressure,  a  pressure  drop  leads  to  a  longer  re¬ 
sponse  time  than  the  corresponding  pressure  rise  -  because  the  equilibrium 
process  takes  a  longer  time  near  the  lower  final  pressure.  Thus,  the  plots 
presenting  phase  lag  show  two  points  per  run  per  frequency  as  an  indication 
of  the  range  of  variation.  Both  points  were  obtained  at  the  piston  midstroke 
position  (or  zero  normalized  pressure  in  Figure  42)  --  one  180  degrees  from 
the  other. 

Figure  43  presents  a  sample  readout  from  the  flush-mounted  piston  trans¬ 
ducer  and  illustrates  the  fact  that  the  pressure  did  not  vary  sinusoidally. 

By  assuming  Boyles's  Law  to  be  valid  here,  it  can  be  shown  that  the  non- 
sinusoidal  variation  is  to  be  expected: 

PV  =  P0  V0 
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where  p  =  pressure 
V  =  volume 

subcript  0  denotes  condition  at  piston  miJstroke 

The  volume  at  any  instant  can  be  expressed  in  terms  of  the  angular  position 
of  the  crank  0, 


V  =  V0  +  AV  sin  0 
=  V0  (1  +  V1  sin  0) 


where  V'  =AV/V0 

Substituting  this  value  of  V  into  the  original  equation, 

p  =  p0/  (1  +  V'  sin  0) 

Expanding  and  neglecting  higher  order  terms  of  V',  the  expression  can  be 
written 


P  "  P0  =  "  Po  V'  sin  0(1  “V  sin  0) 

Note  that  as  V'  becomes  small,  p  -  p0  approaches  -  p0  V  sin  0.  However, 
for  this  investigation, V'  was  not  small. 

QO 

A  report  of  some  interest  is  that  of  Bergh  and  Tijdeman  ,  wherein  are 
presented  results  of  theoretical  and  experimental  investigations  of  the 
relationship  between  the  sinusoidal  pressure  disturbance  in  a  given  volume 
and  an  adjoining  volume.  However,  the  theoretical  approach  employed 
very  small  sinusoidal  disturbances.  The  present  investigation  employed 
disturbances  of  an  order  of  magnitude  greater,  more  realistically  simula¬ 
ting  the  disturbances  in  the  wind  tunnel. 

Examining  Figures  38  through  41,  the  following  trends  are  noted: 

1.  The  pressure  attenuation  peaks  (first  resonance  peak)  occur  at 
lower  frequencies  as  the  tube  length  is  increased,  and  the  peak 
values  decrease  with  increasing  tube  length. 

2.  The  phase  lag  increases  with  increasing  tube  length  and  with  in¬ 
creasing  frequency. 

3.  The  amplitude  of  the  oscillating  pressure  has  a  significant  effect 
upon  pressure  attenuation  and  phase  lag  --  the  effect  being  more 
pronounced  with  the  greater  tube  length. 
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In  addition  to  trends  (1)  and  (2)  stated  above,  Bergb  and  Tijdeman  indicate 
that  pressure  attenuation  and  phase  lag  are  functions  of  pressure  level. 
Pressure  level  was  not  varied  during  this  investigation.  Trend  (3),  how¬ 
ever,  is  not  noted  by  them,  probably  due  to  the  small  magnitude  of  pres¬ 
sure  disturbuance  that  they  imposed. 

The  results  of  the  short  (6-inch)  tubing  investigation  are  presented  in 
Figures  44  and  45.  It  is  immediately  obvious  that  connecting  to  the  plug 
side  of  the  transducer  produces  greater  phase  lag  and  greater  deviation  of 
attenuation  from  1.  00,  probably  because  of  the  larger  volume  between  the 
tubing  and  the  diaphragm.  Within  the  frequency  range  (up  to  32  cycles  per 
second),  mean  amplitudes  (61  and  31  inches  of  water),  and  tube  diameters 
(.049  and  .  094  inch)  tested,  phase  lag  and  attenuation  appear  to  be 
nearly  independent  of  amplitude  and  diameter  for  the  flush  diaphragm  face 
connection.  The  pressure  attenuation  varied  from  1.  00  to  1.  02  nearly 
linearly  with  increasing  frequency,  and  the  remote  pressure  lagged  the 
actual  pressure  by  an  amount  that  varied  from  0  to  2-  1/2  degrees  nearly 
linearly  with  increasing  frequency.  Also.the  phase  lag  showed  only  a 
slight  variation  over  a  cycle. 

CONCLUSIONS 


Since  the  scanivalve  gasket  restriction  of  the  ports  has  a  strong  effect 
upon  the  data,  the  use  of  the  scanivalve  and  single  transducer  is  in  all 
likelihood  precluded. 

In  addition  to  being  functions  of  tube  geometry  and  frequency  of  oscillation, 
pressure  attenuation  and  phase  lag  are  also  functions  of  pressure  level, 
amplitude  of  pressure  oscillation,  and  portion  of  cycle.  Therefore,  it 
would  be  a  very  difficult  task  to  determine  the  instantaneous  local  oscilla¬ 
ting  pressure  from  remote  readings. 

For  the  above  reasons,  it  is  highly  unlikely  that  a  scanivalve  and/or  trans¬ 
ducers  mounted  outside  the  wind  tunnel  car.  be  employed.  It  would  be 
necessary  to  mount  transducers  inside  the  model,  where  there  would 
still  be  need  for  short  tubing  to  connect  some  of  the  orifices.  This  connec¬ 
ting  tubing  should  be  kept  as  short  as  possible  and  should  not  exceed  the 
6- inch  length  tested.  •  By  directly  using  the  remotely  read  data  of  the 
short  tubing,  the  ensuing  errors  would  be  within  acceptable  bounds.  The 
°rror  level  could  be  reduced  by  assuming  an  average  or  constant  lag 
ihroughout  the  cycle,  shifting  the  oscillating  pressure  curve  the  appropri¬ 
ate  amount,  and  applying  an  attenuation  factor  (which  would  be  the  recipro¬ 
cal  of  the  attenuation). 
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PRESSURE  ATTENUATION 


Tubing  length 


1  1 

1 

Different  sym 

:  denote  different 

runs  at  the  . 

1  1 

•ie  conditions 

1  1 

FREQUENCY  -  HERTZ 


FIGURE  38.  Pressure  Attenuation,  Influence  of  Tube  Length,  Using 
Scanivalve,  Mean  Pressure  Amplitude  PI  in.  H2O,  Tube 
Diameter  0.  04  9  in. 


PHASE  LAG  -  DEGREES 


0  4  8  12  16  20  24  28  32 

FREQUENCY  -  HERTZ 

FIGURE  39.  Phase  Lag,  Influence  of  Tube  Length,  Using  Scanivalve, 

Mean  Pressure  Amplitude  61  in.  HgO,  Tube  Diameter  0.049 
in. 
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PHASE  LAG  -  DEGREES 


FIGURE  39.  Continued* 
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FIGURE 


16  20  24  28  32 

JENCY  -  HERTZ 

9.  Continued. 
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PRESSURE  ATTENUATION 


FIGURE  40.  Pressure  Attenuation,  Influence  of  Tube  Length  and  Pres¬ 
sure  Amplitude,  Tube  Diameter  0.049  in. 
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FIGURE  41.  Phase  Lag,  Influence  of  Tube  Length  and  Pressure  Ampli¬ 
tude,  Tube  Diameter  0.  049  in. 
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PRESSURE  ATTENUATION 
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FREQUENCY  -  HERTZ 


FIGURE  44.  Pressure  Attenuation,  Short  Tubing  Length  (6  in. ) 
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PHASE  LAG  -  DEGREES 


FIGURE  45.  Phase  Lag,  Short  Tubing  Length  (6  in.). 
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APPENDIX  III 

INSTANTANEOUS  PRESSURE  COEFFICIENTS 


For  the  instantaneous  pressure  coefficient  data  sheets,  the  symbols  are 
defined  as  follows: 

AA  Angle  of  attack,  a 

DELTA  AA  Oscillating  Amplitude,  Aa 

MEAN  AA  Mean  angle  of  attack,  a 


NOTE: 


X/C  equals  .01  is  really  .007  5.  This  was  rounded  off  to 
.01  for  tabulation  purposes  only. 
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INSTANTANEOUS  PRESSURE  COEFFICIENTS 
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INSTANTANEOUS  PRESSURE  COEFFICIENTS 
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INSTANTANEOUS  PRESSURE  COEFFICIENTS 
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RUN  NO  55  K  ■  *1897  DEL T  AA  ■  6*11  NEAN  AA  “  -.20 

.20  1.38  2.85  4.12  5.09  5.70  5.90  5.70  5.09  4.12  2.85  1.38 
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